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INTRODUCCIÓN 


En el presente texto de analizan los 
tipos de transformaciones integrales más 
frecuentemente utilizadas en la solución de 
problemas de la mecánica de la construc- 
ción: transformaciones integrales de Lapla- 
ce, Fourier, Hankel y Mellín, y se presen- 
ta la solución de diferentes problemas está- 
ticos y dinámicos haciendo uso de las 


mismas» 


El libro está dirigido a los estudian- 
tes de las Facultades de Ingeniería Civil e 
Industrial y también a ingenieros de estas 
especialidades que se dedican al trabajo de 


investigación. 


Con el incremento de la complejidad de 
las soluciones constructivas de las obras y 
los edificios, a la cual está vinculada tam- 
bién la exigencia de una mayor seguridad y 
durabilidad de las mismas, se hace cada vez 
más necesario la obtención de las solucio- 
nes exactas de las diversas tareas de la me- 
cánica de la construcción, las cuales, no 
solo poseen interés por sí mismas, sino ade- 
más, porque sirven para la comparación de 
las soluciones alcanzadas por diferentes 


métodos aproximados. 


Uno de los métodos más efectivos de ob- 
tención de soluciones analíticas de diferen- 
tes tareas de la mecánica de la construc- 
ción es la aplicación de las transformacio- 
nes integrales! con su ayuda se pueden re- 
solver las ecuaciones diferenciales a que 
en muchos casos nos llevan los problemas de 
la Mecánica de la Construcción, en deriva- 
das ordinarias o parciales, logrando dismi- 
nuir en estas últimas el número de varia- 


Lles independientes o transformando, en 


amlos casos, las ecuaciones diferenciales 


en ecuaciones algebraicas. 


La posibilidad de la utilización conjun- 
ta de las tablas de las transformaciones 
directas e inversas permite formalizar el 
procedimiento de solución de las ecuaciones, 
y de esta misma forma, aumentar la efecti- 


vidad de su aplicación, 


En la presente publicación se analizan 
las transformaciones integrales más frecuen- 
temente utilizadas, y se presenta además, 
las soluciones de una serie de tareas estíá- 
ticas y dinámicas de la Mecánica de la Cons- 


trucción sobre la base de su utilización, 


En el primer capítulo se presenta la de- 
finición de la transformada de Laplace y se 
describen sus propiedades. Después sobre 
la base de la transformación de Laplace se 
expone un método general de solución de 
ecuaciones diferenciales por medio del cual 
se Obtienen las ecuaciones del método de 
los parámetros de origen para el caso de la 


viga sobre soportes elásticos. 


Se explica la solución de la tarea sobre 
las oscilaciones de un sistema de un solo 
grado de libertad, el problema de las vibra- 
ciones longitudinales de una barra, así como 
de las vibraciones transversales forzadas 
de una viga. Se muestra, además, la posi- 
bilidad de solución de las ecuaciones in- 
tegrales de Volterra, hacia las cuales con- 
ducen las tareas del cálculo de estructuras 


con propiedades reológicas, 


El segundo capítulo se dedica a la apli- 
cación de las transformaciones de Fourier, 
se define y se describen sus propiedades. 
Luego, con ayuda de esta transformación se 
resuelven las tareas de la flexión de vigas 
de longitud infinita y semiinfinita apoya- 
das sobre un medio elástico de Wincler, de 
placas semiinfinitas y franjas semiinfini- 


tas, además de determinarse los valores de 


las fuerzas interiores para el caso del 
estado tensional plano. En el epígrafe 
final de este capítulo se presenta el méto- 
do general de solución de diferentes tareas 
con la ayuda de las transformaciones fini- 


tas de Fourier. 


En el tercer capítulo se introducen las 
transformaciones de Hankell con ayuda de las 
cuales se resuelven las tareas en coordena- 
das polares, que se presentan en los casos 
de análisis de regiones circulares. Con 
ayuda de las transformaciones de Hankell se 
presenta la solución de la tarea sobre la 
distribución de las tensiones en los cuer- 
pos de revolución en los casos de existen- 
cia de simetría axial. Luego se determi- 
nan los modos propios de vibración de una 
placa infinita. Y por Último, con ayuda de 
las transformaciones finitas de Hankell se 
determina la distribución de la temperatu- 


ra en un cilindro circular. 


Para finalizar, en el cuarto capítulo se 
definen y describen las propiedades de la 
transformación de Mellin con ayuda de la 
cual se resuelven las tareas de cálculo de 
las construcciones en regiones en forma de 
cuñas (cuneiforme). En calidad de ejemplo 
de aplicación de estas transformaciones se 
presenta la solución de la tarea de deter- 
minación del estado tensional plano en una 
región cuneiforme y de una placa cuneifor- 


me comprimida, 


NOTA DEL TRADUCTOR 


Al igual que la monografía El método de 
los elementos finitos en las aplicaciones a 


la teorta de las contrucciones del profesor 


N.N. Leontiev y I.I. Diomin; el presente 
trabajo Aplicación de las transformaciones 
integrales a la solución de las tareas de 
la Mécanica de la Construcción de V.I. Tra- 
vuch, aparecido recientemente, está dirigi- 
do a la preparación de los ingenieros civi- 
les que en sus distintas especializaciones 
se forman en el Instituto de Ingeniería de 
la Construcción de Moscú (MISI), y su tra- 
ducción al español pretende contribuir al 
esfuerzo de acercar los mismos al profesor 


y al estudiante cubano de esta especialidad. 


Desde hace ya algunos años, el tema de 
la transformada de Laplace forma parte de 
los cursos de matemática superior que se le 
imparten a los estudiantes de Ingeniería 
en Cuba. 


Sin embargo, el enfoque por una parte 
más general, desde el punto de vista de la 
presentación de varios de los numerosos mé- 
todos de transformaciones integrales apli- 
cadas actualmente en diferentes ramas de 
la matemática y las ciencias aplicadas, 
incluida la transformación de Laplace y 
por otra parte más específico de esta pu- 
blicación, que está dedicada a los estu- 
diantes de las especialidades de la cons- 
trucción con ejemplos especialmente selec- 
cionados para ellos, explican el interés 


que nos motivó a realizar su traducción. 


La monografía, fue realizada también 
como una tarea extraplán en el ámbito de 
las actividades de saludo al 67 Aniversa- 
rio de la Revolución Socialista de Octubre 
y al 28 Aniversario del Desembarco del 
Granma y de la Fundación de las Fuerzas 


Armadas Revolucionarias. 


Ing. Angel Emilio Castañeda Hevia 


Capítulo 1 


TRANSFORMACIÓN DE LAPLACE 


1.1 DÉPINICIÓN DE LA TRANSFORMACIÓN 
DE LAPLACE 


Durante la solución de ecuaciones dife- 
renciales relativas a la funcion f(t) con 
argumento variable desde t =0 hasta tf =o00, 
se utiliza frecuentemente la transforma- 
ción de Laplace, la cual transforma la fun- 
ción dada f(t) en una función F (P) de la 
variable compleja p = 0 + 17 de acuerdo con 


la fórmula: 


Pb). de fire Ptas (1.1) 


donde se supone que la función FUEJ=0 
para (t< 0), y además, que la función es 
continua conjuntamente con sus n Primeras 
derivadas, con la excepción de algunos pun- 
tos aislados en los cuales la derivada 

f" (£) permite discontinuidades de primer 
género. Además, se supone que la función 
f(t) al tender a infinito satisface la con» 
dición de no crecer más rápidamente que una 


cierta función exponencial con Índice de 


crecimiento S, de forma tal que: 


S,t 
IF(é)N <Me (12) 


La condición f(t) = 0 para (t < 0), ge- 
neralmente, se satisface en diferentes pro- 
blemas físicos con condiciones iniciales 
dadas para un valor determinado del tiempo, 
el cual nosotros siempre podremos adoptar 
como origen del cálculo t = 0. La función 


más simple que satisface esta condición es 


la función de Heaviside: (o función paso 


unitario) 
te > 0) 

n(t) = (1.3) 
o (t<o) 


Cualquier función p(t) multiplicada por 
la función n(t) también satisface esta con- 


dición. 


La función F (p), determinada por la 
ecuación (1.1) se llama también imagen de 
la función f(t) o su transformada. A veces, 
con el objetivo de simplificar la escritura 
de la ecuación en lugar de (1.1) se utiliza 


la nomenclatura siguiente: 
F(p) = L[f (£)] 


Si se conoce la función F(p), la origi- 
nal f(t) se determina por la fórmula si- 


guiente: 


a+t 0 


ARA (1.4) 


at aio 
En la ecuación (1.4), la integral se 
toma a lo largo de la recta Re, =a> So, 
la cual garantiza su convergencia y se com- 
Prende como el límite de la integral a tra- 
vés del segmento (a - 1h), (a + Tb) cuando 
b=>>% (hb tiende a infinito). 


Determinemos en calidad de ejemplo la 
transformación de Laplace de la función 
nN(t). Sobre la base de (M1.T) y (1:23): 


1 


er near (1.5) 


Si: 
pp? Py? 
vit)=e ,fit)=e ns) 


entonces: 


oo (P, -P)Jt 1 
Ep) =5 nít)e dt =—— (1.6) 


1.2 FROPIEDADES DE LA TRANSFORMACIÓN pk 
DE LAPLACE po -k 


1.2.1 Propiedades de linealidad 1 
UR) - — E - — A 


De la ecuación (1.1) se comprende que si p -k p + kx 
f, (6) y f, (6) son dos funciones, que tienen 1 
por transformadas F, (p) y tes (p) y 4 y É son = _—Á 
constantes, entonces: E 
af, (6) + Bf, (611 =0F, (9) +BF,() 0.) UVR] == R E — == k => - 
p" -k p" +k 
Esta propiedad permite determinar la 
transformada de una función que está com- sie KE . 
p -Kk 


puesta por la suma de otras funciones cuyas 


transformadas son conocidas; por ej lo 
EP 1.2.2 Propiedad de semejanza (o del cambio 


la ecuación (1.7). 


de escala) 
1 twt —Lt 
L[ sen wt] = Ll 0 (e e JJ = Determinemos la transformada de la fun- 
ción f (at) donde 4 es una constante (4 > 0): 
1 1 1 
21 ( p- tw p +10 pa 50 £ 
flat) =f f(at)e P dt = 
= Se — e T 
2 on 
port $ = E ¿ e S E F L, 
e qa o e E A (1.8) 
De forma semejante tenemos la ecuación 1.2.3 Tr. EÓ 506 taa 
siguiente: 
Determinemos la transformada de Laplace 
Ll cos we] = p o L[ senh «t] = de la función diferencial Fue , donde 7 
p AY) es el orden de la derivada. 
eS sin = 1, integrando por parte obtenemos: 
, p — me po 
pres Pat [fi Pd + 
L[ cosh wt]= A y co -pt 
pp -o +p5 fi) e Pide =-f(0) + FI) 


Si utilizamos las expresiones anterio- . 
En este caso se tuvo en cuenta la condi- 


s demo cribir 1 £ j Y =- d . 
AS A O IS ción (1.2) para la evaluación de límite 


formadas de A.N. Krilov: 
5 pes E cuando Ll >", 


aa Me =p pe d por de De forma semejante realizando integra- 
ciones sucesivas por partes llegamos a la 
E do: AN dependencia buscada: 
4 4 
E $ ES 
are =p Fw - po fo) - peo) 
ara. ——_— + —» $8 
AA AE A TR (1.9) 


donde: 
FO AUD o, 


Corresponden a los valores de la función 
y de sus derivadas para t+ = 0, 
1.2.4 Derivada de la transfonuada 


Derivando de forma inmediata ambas par- 
tes de la igualdad (1.1) obtenemos que: 


- 


p) (p) = en” y AR 
a partir de lo cual: 


4 cat )] =P p) (1.10) 


Si tomamos como ejemplo si f(t) =nm(t) ” 
a partir de la ecuación (1.5); 


1 
£E[ 1] e 


y en este caso: 


na! 


n 
A 
Pp 


po) (p) 


consecuentemente y de forma directa obtene- 
mos la imagen t” 

q os 
1 


1.2.5 'lransformada de la integral 


Determinemos la transformación de la fun- 
ción g(t) donde: 


gt) =5%f (0) de 
de forma tal que: 
g tr = fit) 


lo que aplicando la transformación de La- 
place a ambas partes de esta igualdad obte- 
nemos que: 


PG(p) = F(p) 


a partir de lo cual: 


L[(t)1 =F(p) p* (1, 11) 


o lo que es lo mismo la transformada de la 
integral implica el cociente de la función 


por la imagen. 


1.2.6 Integral de la transformada 


Calculemos: 
S F (p)dp $ PS f(t)e E “de 
o lo que es lo mismo: 


af (t). +7] = [“F(prap (1.12) 


Consecuentemente, la integración de la 
transformada es equivalente a la división 
de la función original sobre su argumento; 
por ejemplo de (1.6) tenemos que: 


bt at 1 1 
in a JS 


a partir de lo cual: 


E A O 
t b GAO 
e ptr 
n= hb 


1.2.7 Teorema del retraso 


Encontremos la imagen de la función 
ft - 71): 


7 - me Pude = Ft - 110 Poe 
Tf, 


cea u)e PY - 


ep (p) 


De esta forma, si la función original 


presenta un retraso en un valor 7 entonces: 


auf - 1) =P E(p) (1.13) 


Encontremos, por ejemplo, la transforma- 
da de una función escalonada como la mostra- 


da en la figura 1.1. 


Figura 1.1 


f(t) = Cln(c) + n(t - 7) + n(t - 27) + ..) 


Utilizando las fórmulas (1.5) y (1.13) 


obtenemos: 
-pT - 2pT 
1 e e 
L t = ll — + —— + —___— $+ o... = 
[£N = AS ; 5 ] 
=z C —— 
p(1 - eP*) 


1.2.8 eorena del desplazaniento 
en el dominio n 
Determinemos la transformada de la fun- 
ción: 
p,* 
e ft) 


Luego de sustituir directamente esta 
función en la ecuación (1.1) no es difícil 
convencerse de que: 


aer Fr it)l = Fi -p,) (1.14) 


o lo que es lo mismo, la multiplicación de 


la función original por Pot implica el des- 


plazamiento de la transformada en un valor 
igual a p,. Este teorema permite a través 
de las transformadas conocidas, obtener la 
transformada de esas mismas funciones, cuan- 
do están multiplicadas por una exponencial, 


por ejemplo: 


1 


y e Pot sen we] = 
(p + p, Y + y 


1.2.9 Teorema sobre la nultiplicación 
de las transformadas 


En primer lugar, introduzcamos el con- 
cepto de convolución de funciones. Se co- 
noce como convolución de dos funciones a la 
integral del producto de dos funciones cal- 


culada por la fórmula: 
$ 
E AN ES 
= 1% (T)f (t - rjar 
0. =2 1 


Ahora encontramos la imagen de la convo- 


lución, o lo que es lo mismo, calculemos: 
ot L £ ps 
S e dt / E, mm, (t - rjar 
O e Pr = 
P4, (ridrf e 3 (t - Tidt 


Desarrollando en el interior del inter- 
valo del término de la derecha el cambio 
de las variables t, = l -T, Obtendremos de 
forma definitiva que: 
-pt 


oo -ptT co “1 A 
J f, (me ar f, (E, Je at, P, WE, (p) 


o lo que es lo mismo, se demuestra el teo- 
rema de E. Borel, el cual expresa que Ja 

imagen de la convolución es el producto de 
las imágenes de las funciones que componen 


dicha convolución. 


4, * f,) = F (PF, p) (1.15) 


1.2.10 Teorema de ¡multiplicación 
de las funciones originales 


A partir de la ecuación (1.1) es fácil 


ver que: 
a e pt 
Af, (t) É, (£)] S É, mE, (t)e dt 


Cambiando en el intervalo de la derecha 
UN (t) por la fórmula de transformación 


(1,4), la expresión anterior toma la forma: 


+7 00 


3 F. (qre dadt = 


a e do 


a+ E E 
SFM (e CDtatida = 


dto 


ul 


2n1 


a+ 


] F, (q)F, (o-qdaq 


— 
20€ aio 


Oo lo que es lo mismo, la transformada del 
producto de dos funciones se determina a 


través de sus transformadas por la fórmula: 


ario 


l E (E, p-pa 


af (107, (01 = ar! 


-,00 


(1.16) 


1.2.11 Teorema sobre el desarrollo 


en serie 


Supongamos que la función original puede 
ser presentada en forma de una serie de po- 


tencias: 


k-1 


(Es 
22 k 
ft) = 21 VAT 


entonces, utilizando la propiedad 4, se ob- 
tiene directamente que: 

C so 

k ¿Ko1 


Cc 
GS ERE 


= k 
A O 


(1,17) 


Las propiedades aquí presentadas de las 
transformaciones de Laplace se utilizan 
frecuentemente durante la solución de pro- 
blemas concretos. En estos casos. no es ne- 
cesario en cada ocasión calcular las trans- 
formadas o realizar la integración en el 
contorno. Las tablas de las funciones ori- 
ginales y de sus imágenes se presentan por 
ejemplo en el libro de V.A, Ditkino y 
A.P. Prudnikov: Transformaciones integrales 
v cálculo operacional, Ed. Nauka, 1974, 


1.3 APLICACIÓN DE LA TRANSFORMACIÓN 
DE LAPLACE O LA SOLUCIÓN 


DE ECUACIONES DIFERENCIALES 


ORDINARIAS 
La transformación de Laplace encuentra 


una amplia aplicación durante la solución 
de ecuaciones diferenciales ordinarias con 
coeficientes constantes, en los sistemas de 
dichos tipos de ecuaciones, de ecuaciones 
en derivadas parciales y en tipos especia- 


les de ecuaciones integrales y otras, 


Tenemos la ecuación diferencial ordina- 


ria con coeficientes constantes: 


n n-1 ps 
a, 4 lt) + Go (6) +... + Q, U (8) = f(t) 


(1.18) 


y la condición inicial 


n-1 
A (0) = 4-1 


u(0) = uy u'(0) =u 3 

Usemos de acuerdo con (1.1) la nomencla- 
tura siguiente U(p) es Ja transformada de 
Laplace de la función u(t) y f(p) es la 
transformada de Laplace de la función f(t). 
Entonces, aplicando a ambas partes de la 
ecuación diferencial (1.18) la transforma- 
da de Laplace y teniendo en cuenta la ecua- 
ción (1.2), obtenemos la ecuación algebrai- 
ca siguiente: 


A(p) U(p) = F(p) + B(p) (1.19) 


donde: 
n 
A(p) =4a,p +4, Pp + 0p +2, 
n n-2 
B(p) = 4, (ap d- 2 ap + + 
ao u la po? + A +. +4,) + 


to tau 

n n-1 
Encontramos a partir de 1.19 la función 

U(p), e invertiéndola después en correspon- 

dencia con la fórmula (1.4) hallamos la 


función buscada: 


ato 
1 S F(p) + B4p) Pap 


u(t) = 
arto qee 


(1.20) 


Para la solución definitiva de la tarea 
es necesario calcular la integral de con- 
torno en la ecuación (1.20). Teniendo en 
cuenta que la integración se debe realizar 
a través del contorno, el cual está com- 
puesto de la recta a - 1%, a +1 yde la 
parte del círculo de radio infinito que se 
encuentra a la izquierda del punto a (ver 
EG. Nerds 


Por el postulado de Jordan la integral 
a través del círculo tiende a cero. Deter- 
minemos, entonces, los valores Pr, que co- 
rresponden a las raíces del polinomio 4 (p); 
si los mismos son enteros, entonces tenemos 
la ecuación siguiente: 
En E PP) + Blp) pj 


k£1 A* > 


20) (1.21) 


y sim, es una fracción de las raíces Pr 


entonces: 
mao 
u(t) = ¿En == 1in a 3 . 
E 


ma E 
- Pp) + BP =p) ed (1.22) 


De esta forma, la aplicación de la 
transformación de Laplace permite de una 
manera suficientemente precisa obtener la 


solución de la ecuación diferencial. 


Con ayuda de la transformación de Lapla- 
ce se puede obtener una fórmula universal 
de la línea elástica de la viga (método de 
los parímetros de origen). Si sobre una 
viga actúan en los puntos 4, las fuerzas 
concentradas Dj Y los momentos Moo y ade- 
más, la misma está cargada por la fuerza 
distribuida q (x) entonces la ecuación de 
la flecha de la viga puede ser obtenida a 


partir de la solución de la ecuación dife- 


rencial: 
" _ _M(x) 
donde: 


M(x) =M, + QT + z P, (e - a,) + 2 4, + 
yA t 


x= 
+ Í q(t) (u - t)dt 


Aquí M, Y Q, son los valores del momen- 
to flector y de la fuerza cortante en el 


extremo izquierdo de la viga x = 0. 


Representemos a través de y(0) y y*(0) 
los valores de la función de la flecha 
y (1) y de su derivada en el punto x = 0. 


Si aplicamos a la ecuación (1.23), la 


transformación de Laplace tendremos: 


x 1 E 
Ur ! A 


(1,24) 


x 
1 = 
Uf altre - srde] = 
LS az 
= Fy(p) ; dm” = E, (P) 


entonces, de acuerdo con (1.9) obtenemos: 


y) = y lo) + yr 0) + 
p 


1 
+ A (p) + ME, (p) + O, 


+ z MEP) + P, (p)] (1.25) 


Ahora será necesario invertir la fun- 
ción Y (p), donde si recordamos que antes 
fue obtenido que: 


PO sd (1.26) 
L(1) ES L(x) ST 


Pp 


consecuentemente, (1) y (1) serán las fun- 
ciones Originales de las expresiones p”* y 
P?. Para la obtención de las funciones 
Originales de las restantes funciones uti- 
lizaremos el teorema de Borrel (propie- 
dad 1.1.9). De la fórmula (1.15) y (1.24) 


se evidencia que: 


x P (P) 
Ulx - u) 1 
E A dl > 
12 
d O) 
ri 
Pp 
E, (p) 
us ña” S “qío) (u - £)dt du] + 
x (í- u) (qu — a,) G¿(P) 
Vil ROA z 
a; El (4) do 
(1.27) 


invirtiendo ahora la expresión (1.25) y 
utilizando las expresiones (1.26 y M2), 
obtenemos la función buscada de la flecha 


de la viga en la forma: 


, Ne (T -u) 
y (5) =y(0) + xy” (0) +45 ET) du + 


zx 
du + % S A du + 


u) 
E 

E an 
in 


GER > 


Para la viga de sección constante, so- 
metida a una fuerza uniformemente distri- 
buida la fórmula (1.28) se simplifica: 


po 
Y (=) = y(0) + xy*(0) + M4, ET + 


2 
te == 
1 
qe 3 0 
ar ¿ee * ET 


1.4 VIBRACIONES DEL SISIEMA CON UN GRADO 

DE LIBERTAD 

Determinemos la ley de movimiento de un 
Sistema de un solo grado de libertad bajo 
la acción de un cierto impulso aplicado de 
forma constante sobre él, de magnitud Q 
aplicado en el momento de tiempo ta = kr 
(K-* 0,1,2 -). La masa del sistema m, es 
una fuerza recuperativa proporcional al 
desplazamiento. La ecuación de las vibra- 


ciones de este sistema tiene la forma: 


dé 


m 


+0t=Q r20 5 (t - kr) (1.29) 


Supongamos que la desviación inicial 
del sistema x(0) = 0 Y la velocidad ini- 
cial x'(0) = 0, 


Aplicando a la ecuación (1.29) la trans- 
formación de Laplace y utilizando la fórmu- 
la (1.9) obtenemos: 

2 kr 
Pop) + y =4)=-Qq, 20 aa 
donde se representa $ = cm”! QQ =Qq*, 
y 2 (p) que es la transformada de Laplace 
de la función x (+), determinada por la 


ecuación (1.1) y teniendo en cuenta la pro- 
piedad delta de las funciones tenemos la 


ecuación siguiente: 


y A de AP 


Resolviendo la ecuación algebraica obte- 


nida, se tiene: 


4 


S  —]kr 
x(p) = a ge 
p 5 Pe k=0 
Q 


(p* +4)(1 se 2) 


y la función buscada x(t) se obtiene a tra- 


vés de la transformación inversa: 


qe? 


dp 
p+rdra e?) 


Determinemos los polos de la función 


bajo el signo de integral. Evidentemente 
ellos serán los puntos: 

: 2kat 
p =0,p.=iw Ls e (k = 1,2y.) 


Si suponemos que T no constituye un nú- 


mero entero, la fracción o enton- 
ces todos los polos indicados son simples 
y de acuerdo con la fórmula (1.21) x(t) se 
determina como la suma de las diferencias 
de las funciones bajo el signo integral: 


x(t) EA cos té - 5) - 


7 


1 
=a cos k ——t) 


2 
kr 


oo 
-.,2 
=1 2 


1 


1.4.1 Vibraciones longitudinales 
de la barra 


Determinemos la ley de variación de los 
desplazamientos y (x,6) de una barra empo- 
trada por un extremo (x = 0) bajo la acción 


de la fuerza armónica Q sen wt, aplicada 
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en el extremo libre de la barra x =l. La 
ecuación de las vibraciones longitudinales 
de una barra de sección normal tiene la 
forma: 


2 
Y yuatp — a Eo TN = 0 
dt? dx? 


(1.30) 


donde a: = EF y? y E, F y k son iguales a 
los correspondientes módulos de elastici- 

dad, área de la sección transversal y masa 
lineal de la barra respectivamente. Si su- 
ponemos en el momento inicial de tiempo 


t = 0, entonces: 


y (5,0) = 0, y¿(5,0) = 0 (1.31) 


Escribamos las condiciones de borde. Te- 
niendo en cuenta que la barra está empotra- 
da en un extremo y(0,t) = O y a partir de 
la ley de Hooke se puede establecer que la 


segunda condición de borde corresponde a: 
y¿(Lyt) =Qr>- ET* sen wt 


Aplicando la transformación de Laplace 
a la ecuación en derivadas parciales (1.30) 
y a las condiciones de borde, obtenemos, 
teniendo en cuenta la fórmula (1.9) y las 
condiciones iniciales (1.31), la ecuación 
diferencial ordinaria en relación con la 


función Y (x,p): 


1:67 00 G 
dx" 


Pp Yi,p) - e (1.32) 


y las condiciones de borde transformadas: 


Y (=p) =0 


dY (x,p) ad QF w 
A 0 


(1.33) 


La solución de la ecuación (1.32) tiene 


la forma siguiente: 


Y (xp) =C, cosh e + C, senh L— y 


Satisfaciendo las condiciones (1.33) en- 


contramos los coeficientes: 


QFaw 
Ep(p? +42) cosh P- 1 


La función buscada y (1,t) se determina 
por la fórmula inversa: 
Y 
QFaw a pe dp 
E(p? + «2 )p cosh qe 


senh 


y(r,t) La +7 


Determinando los polos de la función 


bajo el signo de integral: 
P=0,p=ttw jp, = tral”! (k - 0,51) =t0, 


Obtenemos a través de la férmula (1.21) 
el valor definitivo de la función de los 
desplazamientos longitudinales de la barra 
vibrante. 


1.5 VIBRACIONES FORZADAS DE UNA VIGA 


Analicemos las vibraciones forzadas que 
se producen en una viga, sometida a la ac- 
ción de una carga armónica arbitraria fíx,t), 


sen (wí + 4) con una frecuencia dada w, 


La ecuación de las vibraciones forzadas 


de la viga tiene la forma: 


El yt 2) + m Uy test) SpAl. 
qe ap 


= f(x) sen (we +0) (1.34) 


Representemos ( =m (El)"* y tomemos en 
cuenta que las vibraciones forzadas produci- 


das ocurren de acuerdo con la ley; 


Y(T,t) = p(1) sen (we +) 


donde: 


Y (1): Forma de las vibraciones 


Sustituyendo el valor y (x,t) en la ecua- 
ción (1.34) y representando k* = mu? (El) 
obtenemos la ecuación de las formas de vi- 


bración de la viga. 


IV 4 1 
e (u) -k p() * ET Íf.) 


Para la solución de esta ecuación apli- 
quemos la transformación de Laplace, enton- 
ces, tomando en cuenta la fórmula (1.9) la 


transformación en la ecuación algebraica 
con respecto a la transformada % (p): 


pop) - pero) - pep! (o) - pe”"(0) - p"!(0)- 
á : CA 
-= k b(p) = =T F(p) 
En este caso se representa: 
Pp) = E me Paz y Pp) = 1570) e PE 
De la ecuación obtenida determinamos: 


3 2 
d(p) = p(0) ES ABE p'(0) as e 
p* q ad p 3 ke 


+ p"(0) —L 910) Li 
p A p* - 


PS 119 A 


Elp* - xk) 


Teniendo en cuenta las imágenes de las 
funciones de Krilov: S, T, U, V, y la fór- 
mula (1.15) obtenemos la expresión de la 


forma de las vibraciones: 


, 
P (3) = p(0) Síkc) + L 0-22) + 


mm us 
k? ezr* 


U(kx) + 


+. 2ELO) 
l? 


+ Vik(x - u)]du 
(1.35) 


Si por ejemplo, sobre la barra actúa una 
fuerza concentrada unitaria aplicada en el 


punto x, , entonces: 
fi) =85(x =x,) 


y el último sumando en la fórmula (1.35) 


toma la forma: 


(* ED? Uk( - 2,)) 


1.6 ECUACIONES INTEGRALES DE VOLTERRA 


En la mayoría de los casos de tareas de 
la mecánica de la construcción, teniendo en 
cuenta el tipo de materiales que se anali- 
zan, frecuentemente se llega a ecuaciones 
integrales de Volterra de segundo género: 


t 


ple) + S K(t -=7) y (T)dT = f(t) (1.36) 


donde: 
y (t): Función buscada. 
Kk(t - 71): Núcleo de la ecuación. 


f(t): Función conocida. 


Para la solución de tales ecuaciones re- 
sulta sumamente efectivo la aplicación de 
las transformaciones de Laplace. Represen- 
temos a través de Y (p), k(p) y P(p) las 


transformadas de las funciones y(t), k(t) y 


Je 
ción a la integral en la ecuación (1.36) 


Aplicando el teorema de la convolu- 


llevaremos esta ecuación a la forma: 
Y (p) + K(p) Y (p) = P(p) 


Determinando a partir de aquí la función 
Y(pD) y aplicando la transformación inversa 
de Laplace a la expresión obtenida, encon- 


tramos la función buscada: 


ato 
1 d(p) P? 
y (x,t) O S ; 3 R dp 
TL jo FE (1.37) 


Si son conocidas las raíces de la expre- 
sión 1 +Xk(p), la integral en la ecuación 
(1.37) puede ser calculada en forma clara 
a través de las fórmulas (1.21) y (1.22). 
De esta misma forma puede ser resuelta la 
ecuación de Volterra de primer género. 


t 
Y K(£ -T) y (rjar = f(£) 


En este caso la ecuación buscada toma 


la forma: 


a+to 
E. 1 P(p) t 
pu) O 


De forma semejante se resuelven los sis- 
temas de ecuaciones integrales de Volterra 
o los sistemas de ecuaciones integro-dife- 


renciales. 


Capítulo 2 


TRANSFORMACIÓN DE FOURIER 


2.1 DEFINICIÓN Dri LA TRANSFURMACIÓN 
DE FOURIER 


Durante la solución de diferentes ecua- 
ciones diferenciales en coordenadas rectan- 
gulares se pueden aplicar en forma efecti- 
va las transformaciones integrales de Fou- 
rier, las cuales transforman una región 
continua dada o una función seccionalmente 
continua [f(%:)] que puede ser integrada to- 
talmente en el intervalo (-",%) en una fun- 
ción F(A) de acuerdo con la fórmula: 


A O e (2.1) 


zz 2 


F(A) = 


donde: 
A: Parámetro de transformación. 


Las funciones F(A) se llaman funciones 
transformadas o transformada de Fourier de 
la función f(x), y a veces se representa 
como F[f (2:)]. 


Si es conocida la función F(A), sa fun- 
ción f (wc) puede ser determinada por la fór- 
mula de la inversión (o fórmula de la fun- 
ción inversa). 

1 pe PA) ¿A A 


— 


ME AE 


(2.2) 


f(x) = 


Para la función f(x), simétrica con 
respecto al origen de coordenadas, tenemos 
a partir de (2.1) y (2.2) que: 


FM) = VÍ ST) cos Arde (2.3) 
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fu) = VE ¡7 EL) cos de (2.4) 


Si la función f(x) es inversamente simé- 
trica se Obtiene: 


FA) = a f (<) sen hwdx (2.5) 


2 on 
f(x) = rl F¿(A) sen Mudz (2.6) 
Los pares de fórmulas (2.3) y (2.4), 
(2.5) y (2.6) determinan las correspondien- 
tes transformaciones seno y coseno de Fou- 


rier. 

De forma semejante se introducen las 
transformaciones múltiples de Fourier, por 
ejemplo, en el caso bidimensional, la trans- 
formación de Fourier de la función f'(x,y) 
se llama función: 


FE) =FIf(yl = 


=—í fíewy)ye ASHEY) yy 


pl 


La fórmula de inversión, en este caso, 
tiene la forma: 
e ZÚ( 
Puy) - ram DA e, 
«LO 


En el caso de funciones pares f(x,y) con 


relación a ambos ejes. 


MS = =P f(c,y) cos Ax cos ty dxdy 
Fioy) = ESF Y, ME) cos Ax cost yd 
Para las funciones impares (inversamen- 


te simétricas) con relación a los dos ejes: 


F¿MsE) = ES fun sen Ax sen ¿ y dxdy 
fem) == SÍ PLQ,E) sen Me sent Y dt 


2.2 PROPIEDADES DE LA TRANSFORMACIÓN 
DE FOURIER 


2.2.1 Propiedad de linealidad 


Si f, (2) y f, (1) son dos funciones que 
satisfacen las condiciones indicadas ante- 
riormente y 4% y É son dos números arbitra- 
rios, a partir de la fórmula (2.1) es evi- 
dente de la propiedad de linealidad de la 


transformación, lo que quiere decir que: 


Flaf, (0) + Bf, (0)] =aF[f, (<)] + BF[f, (0)] 


(2.7) 


2.2.2 Propiedad de semejanza 


Si tenemos la función dada f (ax), a par- 


tir de (2.1) es evidente que: 


Af. LL rt, (2.8) 


2.2.3 Derivación de la función original 


Determinemos la transformación de Fou- 
rier de la función Pe) , donde n es el 
orden de la derivada. Analicemos primera- 
mente el caso n = 1 y supongamos que la 
función f(x) es continua. Integrando por 


partes, obtenemos: 


so A A co 
1 Sh bla Lo Sh pue £ 
a y 2 
TA 


Pm a ARO 25 
2 DO 


Aquí se supuso, además, que el término 
fuera de la integral se anula dentro de los 
límites señalados. Generalmente, en las 


tareas prácticas esta suposición siempre 


se realiza. La función f(x) no siempre pue- 


de ser continua, y se pueden admitir dis- 
continuidades de primer género, pero en 
este caso, ellas deben ser tenidas en cuen- 


ta en la evaluación de la fórmula (2.9). 


Supongamos que la función f(x) tiene una 
discontinuidad en el punto x= = 0, y en ese 


caso, en el lugar de la fórmula (2.9)ten- 


dremos: 
1 OS thx 1 
S f'(me a ——— [f(-0) - 
Y 21 PS y 2 


- fF(+0)] A F(A) (2.10) 


Aquí, a través de f(-0) y (+0) se repre- 
sentan los valores de la función f(x) a la 
izquierda y a la derecha del punto x = 0. 
De forma similar pueden ser tenidas en cuen- 
ta las discontinuidades en otros puntos del 


eje (x). 


Extendiendo la fórmula obtenida a una 
derivada arbitraria de orden (n) obtenemos; 
i 
1 oo io 
— / f tale * “de = 
yn 


An LOA o)  Plor7 + 


+ GO”) (2.10 


Si la función f(x) es continua, el miem- 
bro fuera de la integral se anula y se ob- 


tiene que: 
O e (FO 212) 
DO 


Ahora pasemos a la deducción de las fór- 
mulas para la función Pla) en el caso de 
las transformaciones seno y coseno. De 


esta forma, integrando por partes, obte- 


nemos: 
IZ 170 cos Áxdx =- IZ 510) + 
+A OS) (2.13) 
IZ paga) cos Áxdx = - Y 510) = 
(2.14) 


== EH 
e ER 


VÁ JT ""60) cos dude =- Y fo) + 
+ Vo) - MEF) (2.15) 
Y 77" (a) cos dada = Y 2 910) + 


+1 WE 10) + MM E.) (2.16) 


Aquí F (0), o) son los valores de la 
función y de sus derivadas para el caso 


x= 0, 


De forma similar: 


Y 1%" (0) sen Axdx = —A F (A) (2.17) 
0 E 


v 1770) sen Axrdx V—Ñ 40) - 


'I 


2 59) 
- ELA) (2. 18) 
a sen Axdw = YH a5* 10) + 
+ PE, A) (2. 19) 


Ia sen Made = YH 17"10) - 
2) 3 
ss Y ÑA f(0) + MEA) (2.20) 


Si la función dada f(x) es continua desde 
(+ <x<o), los miembros fuera de la in- 
tegral en las expresiones (2.13) a (2.20) 
se anulan, 


Sin embargo, durante la resolución de 
tareas para las regiones, las cuales se 
determinan en el espacio semiinfinito 
(0<x<o*%), fuera de las integrales en 
estos casos, es necesario tener en cuenta, 
estos miembros debido a que en este caso 
ellos constituyen valores de borde la fun- 
ción f (1) y de sus derivadas. 


De esa forma, y a partir de las fórmu- 


las anteriores, resulta evidente que en el 


caso de estar dadas en el borde los valo- 
res de las derivadas pares, por ejemplo 
f(0) y f"(0) es necesario utilizar en la 
solución de la tarea la transformación- 


seno. 


Si en la misma frontera estuvieran 
dados, los valores de las derivadas impa- 
res f'(0), f”"'(0) se aplicarían las trans- 
formaciones coseno. En estos casos, se 
extiende-la función f(x) dada en el inter- 
valo semiinfinito a través de todo el eje 
con la introducción de derivadas generali- 


zadas. 


En la teoría de las funciones generali- 
Zadas se demuestra que si la función f(x) 
tiene en cualquier punto dado una discon- 
tinuidad de primer género, su derivada ge- 
neralizada será igual a la suma de la de- 
rivada ordinaria (esto quiere decir, la de- 
rivada obtenida sobre la base de la supo- 
sición de la continuidad de la función f(x) 
y a la magnitud de la discontinuidad multi- 
Plicada por la función delta, esto es 
que: 


Fu = f(x) + 5, (2 - 2) 


Aquí f'(x) y f'(<) son las derivadas 
generalizadas y ordinarias de la función 
(1), S, es la magnitud del salto de la 
Y 2 


función en el punto de discontinuidad X,. 


Esta fórmula puede ser extendida para 
el caso de derivadas de orden superior to- 
mando en cuenta los valores correspondien- 
tes de los saltos en el punto de discon- 
tinuidad S¿* 

EZ 


1 p 
RRA 
¿Ey 8/5 (a - xy) 


F(k) 


PO y PO pc) 4 


(2.21) 


A las funciones f(x) y a sus derivadas 
F (K) (1), extendidas ahora a través de toda 
la longitud del eje se pueden aplicar las 


transformaciones de acuerdo con las fórmu- 
las ka TN (2.3) o (2.0) 


2.2.4 'eorema sobre el deslizamiento 


Encontremos la transformación de Fourier 
de la función (% - x,) donde x, también 


puede ser variable tenemos: 


1 00 Az, 
——f fi. -x2,)e da 

Zn pa] 

tx : 
a. e ; S fia “de 
2 DO 
De esta forma: 
TAx, 

F[f(r= -x,)] =e FA) (2.22) 


Si es conocida la transformada F(A) de 
la función f(x), entonces, la transformada 
de 


la 


la función f(x - x,) se determina por 
fórmula (2.22). 


2.2.5 eorema ue la multiplicación 
de las transformadas 


Sean F, (A) y E, (A) las transformadas de 
las funciones f, (u) y Fs (0). Encontremos 
la transformación de Fourier de su produc- 
to, entonces; 


A GRO Di 
LE a 


1 Mu a, 


1 0 O 
- A E, AjJe" "dk PR (u)e 


pan 


A 


an] 


A e 
A y Sa 03 


A (u)f, (0 - urdu 


vV 2n PU 


Lo que da por resultado las funciones 


siguientes: 


——s f, ()f, (x — u)du = 
Va 
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1 poo 
= —— j (u) (x - u)du 
Ta ff, 


Se llaman convolución de las funciones 
f, (u) y f, (1), los cuales se representan 
también como f, (%) Y f, (0). De esta foxma, 
se demuestra el teorema de E. Borrell para 


la transformada de Fourier, 


FF), 2 01=f(0* 7) (2.23) 
1 2 1 2 


Esta fórmula permite calcular, por ejem- 
plo, las integrales de los productos de las 
transformaciones a través de la convolución 


o a la inversa. 


2.2.6 Igualdad de Parseval . 


Analicemos: 


SF, (A) E, Qda = 


DO 


=P MAS E as 
0 


27 LO 
= 17 f, (u)du or ae du = 
an > o ? 
=j” f, (u) f, ujdu (2.24) 
DO 


N 

Esta igualdad puede ser útil durante el 
cálculo de integrales del producto de las 
transformaciones de dos funciones, o por el 
contrario, en el caso en que sea necesario 
calcular la integral del producto de fun- 
ciones f, (u) y f, (-u) puede ser posible en- 
contrar esta igualdad solamente a través 
de la integral de los productos de las 


transformadas. 


En el caso de funciones pares la igual- 
dad de Parseval toma la forma: 
E si £, (u)d 
[Fs 0) FOO =[ f,(u) f, (uydu 


(2.25) 


Y las impares: 


de E A) F,,AJdA = 5 Ff, (uf, (uddu (2.26) 


2.3 FLEXIÓN DE UNA VIGA INFINITA SOBRE 
UNA BASE ELÁSTICA 


En calidad de primer ejemplo de aplica- 
ción de las transformaciones de Fourier 
analicemos el caso de la flexión de una 
viga suficientemente elástica que descansa 
sobre una base elástica en un medio de 
Wincler y sometida a una carga arbitra- 


ria q (uy). 


Si la carga se encuentra aplicada a una 
gran distancia de los bordes, entonces du- 
rante la determinación de las deformaciones 
y de los esfuerzos interiores en la viga en 
regiones cercanas a la Carga se puede des- 
preciar la influencia de sus bordes y ana- 
lizarla a través del esquema de cálculo de 
una viga de longitud infinita, 

La ecuación de la flexión de la viga es 


la siguiente: 


dy, (kr, ) 
das 


(2.27) 


+Ky"(2,) =q,(x,) 


Puede ser transformada en una ecuación 
de coordenadas adimensionales x = TA 
(donde a = (k/4zI)'“, y = k,B, 


donde: 


ko: 


(7 
B: 


Coeficiente de pastel. 


ancho de la viga. 


El: Rigidez., 


Y con las nuevas funciones y, (x, a) , 
q (5) = q, (,a) (El)? , obtenemos la ecua- 


ción: 


(2.28) 


Ese. 4 y (2) = q(e) 


Para la solución de la ecuación obtenida 
aplicamos a ella la transformación de Fou- 
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rier, representando: 


Ax 


A y (u)e* 


YA) = dee (2.29) 
v 27 DO 
QA) = ama (2.30) 


1 
vV 2 pa 
Entonces, teniendo en cuenta la propie- 
dad (2.12) encontraremos a partir de (2.28) 
que: 


qa) 
A EA 


Y (A) = (2.31) 


Transformando esta expresión por la fór- 
mula (2.2) determinamos la flecha de la 


viga: 
SS —lAx 
A A CN (2.32) 
VA + 1 


Analicemos el caso particular de la 
carga en forma de fuerza concentrada p, 
aplicada en el centro de la viga: 


q, (5,) =0p8 (x,) 


Durante el cambio a coordenadas adimen- 
sionales y a las nuevas funciones obtene- 
mos que: 


q(a) = 4 pa 5 (x) 


Consecuentemente, de la fórmula (2.30) 


se obtiene: 


4 pa 


qu) = 
Y ak 


Sustituyendo esta expresión en (2.32) ob- 
tenemos: 


“e cosAx 
nas 


y) - LE y A 


Calculando la integral, encontramos la 
fórmula buscada para la determinación de la 


flecha de una viga de longitud infinita so- 


metida a la acción de una fuerza concentra- 


da aplicada a ella: 


Pare 


a 
e “(cos + sen 1) 


En el caso analizado de carga, teniendo 
en cuenta la simetría del problema, la solu" 
ción de la tarea pudo ser obtenida directa- 
mente aplicando a la ecuación (2.28) una 
transformación coseno. Si en el centro de 
la viga se aplica una fuerza en forma de 
momento concentrada M, entonces: 


q (5,) =M0'(xmÍ,)+. 27 (2.33) 


Y 


qíu) =07M5'(x) . (2k)” 


Esta carga es antisimétrica relativa al 
origen de coordenadas y para la solución de 
la ecuación (2.28) será necesaria la apli- 
cación de una transformación seno de Fou- 


rier. Representemos en este caso: 


Y A) = VHS y) sen ÁArdx 


l (2.34) 


Dd 


QM) a sen Medx (2.35) 


Obtenemos de la ecuación (2.28) la si- 


guiente: 


qm) 


Y (A) = 
2 NETA 


Transformando esta expresión por la fór- 


mula (2.26) y tomando en cuenta que: 


QA) 


2 AN 
VALDE (o) sen ddr 


AM 


y 2k 


encontramos la expresión para la flecha de 


la viga bajo la acción del momento exterior 
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aplicado a ella (M): 


2 oo 
y (o) = - HE ¡%Asenmde y 

Ú A +4 
Si calculamos la integral, obtenemos: 


2 —- 
y (5) = - — 55 


sen zx (a > 0) 
2.4 FLEXIÓN DE UNA VIGA SEMIINFINITA 
SOBRE UNA BASE ELÁSTICA 


Si la viga, que descansa sobre una base 
elástica se encuentra cargada cerca de uno 
de sus extremos, entonces durante la deter- 
minación de sus deformaciones y sus esfuer- 
zos interiores se puede aplicar el esquema 


de cálculo de una viga semiinfinita. 


Para la solución de la ecuación (2.28) 
en la cual 0 <x <* en dependencia de las 
condiciones de borde dadas en la frontera 
se pueden aplicar transformaciones-seno o 
transformaciones-coseno de Fourier. De 
esta forma si las condiciones de borde en 
el extremo corresponden a un apoyo articu- 


lado, entonces: 


y" (0) = y (0) =0 (2.36) 


En este caso, como se deduce del análi- 
sis de la propiedad (3.1), para la solución 
de la tarea se puede aplicar directamente 
una transformación-seno. Entonces, tenien- 
do en cuenta la fórmula (2.20) y las condi- 
ciones de borde (2.36), encontraremos a 
partir de (2,28) la función y(A) y luego 
por la fórmula de inversión el valor bus- 
cado de la flecha, 


y(x) = Ej 


q, a) sen Ax 
dA (2.37) 


Moa 
El valor AS) se determina por la fór- 
mula (2.5), 


Por ejemplo, sia la viga se le aplica 


en el punto x = 4 una fuerza concentrada PD, 


entonces: 
QM) = VE E [7 (2 — a) sen dr = 


y EL sen la 


Sustituyendo esta expresión en la f6r- 
mula (2.37) encontramos que la flecha de 


la viga en este Caso se determina por la 


expresión: 
= e “[ (cos a + 
+ sen a) cos x senh x + 
+ (sen a - 
- COS a) (sen x cosh 2)] (x< a) 
y (r) = 
eS e” cos x + 


+ sen x) cos a senh a + 


+(sen x-cos )sen a cosh a] (> a) 


De forma semejante, si en el extremo de 


la viga fueron dadas las funciones: 
y'(0) = y'"(0) = 0 


entonces, aplicando a la ecuación (2.28) 
una transformación-coseno y teniendo en 


cuenta la ecuación (2.16) encontramos que: 


ya) 1 Ej 


e (A) cos Ax 
4 


A 


A +4 


En este caso, bajo la acción de la fuer- 


za en el punto a: 
q. = VE cos da 


y la flecna de la viga indica; 


Y) = Lo (cos a + 
+ sen 4) cos x coshx + 


+ (sen q - cos 4)sen < senh 2] (% <a) 
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y(x) = Le (cos x + sen 2)cos A cosh a + 


+ (sen 1 -— cos )sen a senh a] (x > a) 


Como ya señalamos, la solución presenta- 
da sobre la flexión de una viga semiinfinita 
puede ser obtenida solamente para los dos 
casos particulares analizados de condicio- 
nes de borde. 


Sin embargo, utilizando la teoría de las 
funciones generalizadas se puede construir 
un método general de solución de la tarea. 
Con este fin, utilicemos la fórmula (2.21) 
colocando en ella x, = 0yXe4. 
de la sustitución de la expresión obtenida 
en la fórmula (2.28) ella toma la siguien- 
te forma: 


Después 


7 (o) = q) (0) 5,5") — 5,52) - 
- 5,5 '(x) - 5,5 (x). 


Los números 5, representan los paráme- 
tros de origen (esto quiere decir, los va- 
lores de las deformaciones y las fuerzas 


interiores en el extremo de la viga (x = 0): 


donde: 
S, : Valor de la flecha. 
Sy: Ángulo de giro. 


S, » Momento flector. 
S, 2 Fuerza cortante, 


Analicemos por ejemplo una viga con el 
borde libre. En este caso, las condiciones 
de borde toman la forma: 


y” (0) e y'"(0) = 0 


entonces: 


S¿ =S, =0 y la ecuación de la flecha 


de la viga se obtiene como resultado: 


y (0) = q) (2) - 5,5" (0) - 8,8 (2) 


Aplicando para la solución de esta ecua- 
ción una transformacion-coseno y teniendo 
en Cuenta que: 


q, 0) = 17 q (0) cos Mode 


y también: 


ro) =YHE57 y(0) cos kede 


obtenemos: 


Q (1) -S8Q (x) -SQ (x) 

YM) = 0 2 1 3.2 
M +4 

donde: 


q, m) = IES 5 (0) cos Aede =-0 


qe) = 1H 5% 5(0) cos dodo = Y 


De acuerdo con la fórmula de rotación 
(2.4) encontraremos la expresión para la 


flecha de la viga: 


co Y (xr) cos Ae 
2 di + 


2 
Ya el AER 


> 25, po a Z_ di 
A +4 
Como se ve de la fórmula obtenida, el 
primer sumando en ella representa la flecha 
en la viga de longitud infinita, producida 
por la carga q 0), y el segundo sumando 


una corrección; de esta forma: 


y (0) = Y, (1) + s,e * (cos x + sen x)/4 


La constante S, puede determinarse de la 
condición de igualar a cero el momento en 


la viga en el extremos xr = 0, 


M(x) =M, (1) + e “(sen x - cos 2) 


3 

2 

a partir de donde: 
Ss, = 2 M,(0) 


De esta forma, la flecha en la viga de 
longitud semiinfinita: 


y (2) = y, (1) + — UL loJe* (cos x + sen x) 


el momento flector: 
-x 
M(x) =M (1) + M_(0)e "(sen I - Cos 1) 


2.5 FIEXIÓN DE UNA IOSA SEMIINFINITA 


Analicemos la flexión de una losa gran- 
de, rectangular y flexible, empotrada en 
uno de sus bordes (fig, 2.1). Durante el 
estudio del estado tensional y deformacio- 
nal de la losa en las cercanías de la 
carga, aplicada en un lugar suficientemente 
alejado de los bordes libres, se puede des- 
preciar la influencia de estos y analizar 
el esquema de cálculo de la losa semiinfi- 
nita, sometida a la acción de la carga dada 
q (1,4). Como es conocido, las flechas de 
la losa se calcula de la solución de la 


ecuación diferencial en derivadas par- 


ciales: 


4 4 4 
Datolac 4) E 0w pe 0% w SS E 
dat dx* dy? dy 
- o —Á (2.38) 
donde: Ñ 


D: Rigidez cilíndrica de la losa. 


A: Operador de Laplace. 


Figura 2.1 


En el caso dado en el extremo de la 
losa x = 0 se deben satisfacer las condi- 


ciones de borde. 


w(0,y) = w'(0,y) = 0 (2.39) 


Para la solución de la ecuación (2.38) 
aplicaremos a ella la transformación de 
Fourier a través de la coordenada y (* < y 
<>), 


Si representamos: 


w(r,A) = A wa ye ay 
zamr. 


y suponiendo que cuando | yl >“ las funcio- 
nes W, w?, w”, w'" tienden a cero obtenemos, 


teniendo en cuenta la fórmula (2.12) la 
ecuación diferencial: 


A 2 
e 22" +A%) w(x,A) - LED 
de da 
(2.40) 
donde: 
oo LAY, 
Qua) = 5% qye Ydy 
2 bO 


De forma semejante, deben ser transforma- 
das las condiciones de borde (2.39). 
resultado obtenemos que la solución de la 


Como 


ecuación (2.38) debe satisfacer las condi- 


ciones: 


w(0,A) =w'(0,A) = 0 (2.41) 


y también: 


lim (0) 


T> en (x,A) =0;5 


(n = 0,1,2,3) 


(2.42) 


De esta forma, la aplicación de la trans- 
formación integral de Fourier permitió 
transformar la ecuación diferencial en de- 
rtivadas parciales (2.38) en la ecuación di- 
ferencial ordinaria (2.40) con las condi- 
ciones de borde (2.41). 
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La solución de la ecuación (2.40) consta 
de la solución general de la ecuación homo- 
génea y de la solución particular w(x,A), 
esto es: 


-Ax 


A E A E ES 


+ 0(u,A) (2.43) 


como es conocido, la solución particular de 
la ecuación (2.40) debe satisfacer a la 
ecuación homogénea, esto es: 


¿A 


(0,A) =0 5 (tn = Otis 23) (2.44) 


Para su obtención, apliquemos a la ecua- 
ción (2.40) la transformación de Laplace, 
con este objetivo introduzcamos la función: 
- oo -px 
w(p A) Sl wíx,A)e de (2.45) 
y teniendo en cuenta que bajo las condicio- 


nes (2.42) y (2.44) se obtiene la siquiente 


ecuación: 

a A 

A O PA 
z E Pp , 


Obtenemos en el lugar de la ecuación 


(2.40) la ecuación algebraica: 


e? -1P pm - dm 
donde: 


AA) = 17 0, me Pas (2.46) 


De esta forma: 


- 


w(p,A) 5 me 0 se E 
D(p* deL A E 


y la solución particular buscada de la 
ecuación (2.40) se obtiene con ayuda de la 


transformada inversa de Laplace esto es: 


. Xx 
1 ¡an ANA do 


EP qa pp - y 


eii de (2.47) 


Las siguientes aplicaciones las presen- 
taremos por comodidad, sobre el análisis 
de losas concretas cargadas. 

Si a una losa se aplica una carga con- 
centrada en el punto con coordenadas (b,0) 
esto es, 9(2,Y) =DP5 (a - DW (y), entonces: 


pa 


QA) = PY y 5 we ay = 
y 2n 


= —E— 5 (0 - dr 
2r 


Gr) —P [O 8 (e - bre Par = 
Var o 


O  (0<x<bDb) 
pi? par <o) 
y 2D 


Sustituyendo esta expresión en la fórmu- 
la (2.47) obtenemos: 


to ¿Ple - b) 
amm) - —L2 ML 
mi 2n D ae (p - 1)? 


(2.48) 


Utilizando las tablas de las transfor- 


maciones inversas de Laplace, podemos es- 


cribir que: 
0 =(0<x< bh) 
ob, = 
E KQ,x-=b) = (b<x <>) 
2V2DA? 


(2.49) 


En este caso a través de k(A,t) se pre- 


senta la función: 
kA,t) =4A t cosh £ - senh At 
Consecuentemente: 


wi A) = UM + Baxye + 


+ ica maje to zx <b) 
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wit,A) = (A + Haje* + 


+ (CsDhacje 0 y PEO) o <<) 
2V2r DA? 


(2.50) 


Ahora, pasemos a la determinación de las 
funciones desconocidas hasta este momento 
AQ),BM),C(A) yDQA) 


Primeramente encontremos, que: 


lím X(A,t) = —a: - a? 


tí > 

y teniendo en cuenta que I +, w(x,A) > 0, 
entonces, considerando el valor de la fun- 
ción w(x,A), obtenemos que en el infinito: 


AE IS MM A 
42 DA? 


a partir de aquí: 


Qa + e ES 
AQ RL dk y BA) BE 
az Da? aVADN 


Teniendo en cuenta las condiciones de 


borde para x = 0 (2,41), encontremos que: 
CA) = - AA); DA) = - [24(4) + 8()] 


Determinando de esta forma las funciones 
A,B,C y D y sustituyendo ellas en la 
ecuación (2.50) obtenemos el valor de la 
función W(x<,A), y después aplicándole a Él 
una transformación inversa de acuerdo con 
la fórmula (2.2), encontraremos la función 


definitiva de las flechas en la losa w(1,Y)-+ 


2.6 FLEXIÓN DE FRANJAS DE LONGITUD INFINITA 


Ahora analizaremos la flexión de una 
losa lo suficientemente larga para que el 
esquema de cálculo de la misma se pueda 
presentar en forma de una franja infinita 


(figura 2.2). 


Analicemos el caso del apoyo articulado 


de la franja en los lados y = ta. 

En este caso, las flechas de la franja 
se determinan a partir de la solución de 
la ecuación (2.38) con las condiciones de 
korde: 


wx, ta) =0; 4, Es ta) =0 


como es conocido, estas dos condiciones 
son equivalentes a las siguientes: 
dv 


dy? 


ww =0 y =0 


y =ta y = ta 


Supongamos que la franja está cargada 
por una carga simétrica con relación a los 
ejes 1 y Y. En este caso la función de 
las flechas Ww(x,y4) también será simétrica, 
y para la solución de la ecuación (2.38) 
es posible aplicarle a la misma una trans- 
formación-coseno de Fourier a través de la 


coordenada y. Si expresamos: 


(A) =y HH Sota my) cos Ay dy 


Entonces, en el lugar de (2.38) obtene- 
mos la ecuación: 


A E Q4 EA) 
e ) (o, (e/A) A E 


donde: 


QA =/H]17 q (e) cos Aydy 


Transformando también las condiciones ; 


de borde estas toman la forma: 


wíta,A) =0, W'(ta,A) = 0 (2791) 


Tomando en cuenta la simetría de la fun- 
ción de la flecha, la solución de la ecua- 
ción (2.38) se puede escribir en forma de 


suma: 
w(x,A) = A cosh Ax + Bhx senh Ax + 0(x,A) 
(2.52) 
Aquí w(x,A) representa la solución par- 
ticular de la ecuación (2.38) la cual pue- 


de ser obtenida de la misma forma que en 


el epígrafe anterior. 


Por ejemplo, bajo la acción de una fuer- 
za concentrada aplicada en el centro de la 


franja, entonces: 


q(-,y) = pó (x) 6 (y) 


obtenemos a partir de (2,49) que: 


2 k (A,x) 


wlx,A) = 
2 Y 2D 


Satisfaciendo las condiciones transfor- 
madas (2.51) obtenemos el sistema de dos 
ecuaciones para la determinación de las 


coeficientes A y B: 


A cosh Ma + Bad senh a = v(c,A) 


A cosh Ma + B[ 2 cosh ha + ad senh da] = 


= - 1 ?*wa,). 


Utilizando la expresión (2.52) encontra- 


remos que: 


Le pl Md senh Ma cosh ka - Za] 


2 Y 2D MX cost” ha 


(2.53) 


P senh 2 ka 
2 Y 2rD MA cost? ha 


Bx=- 


Si los lados de la franja y = ta están 
empotrados, entonces, satisficiendo las 
condiciones de bordes encontramos el si- 
guiente sistema de ecuaciones para la de- 


terminación de los coeficientes A y B. 
A cosh M + Bok senh da = v(a,A) 


A senh la + B[ senh Aa + adA cosh da] = 


= Au! (a,A) 


La solución de este sistema da: 
mien Pl senh_2%a - A*a*] 
2 Y 2rD MN (senh 2da - 2 ha) 
(2.54) 


-- P senh 2k%a 
2 Y 2D M (senh 2%a - 2 Ma) 


B= 


Las funciones A(A) y £2(A) de las ecua- 
ciones (2.53) o (2.54) sustiuidas en la fór- 
mula (2.52) dan la posibilidad de encontrar 
el valor transformado de la función de la 
flecha en la franja fa, A) para los dos 
tipos de condiciones de borde analizados an- 
teriormente. El valor real de la flecha 
w(2,Y) se determina por la sustitución de 
la función encontrada w(xc,A) en la f£ó8rmu- 
la (2.4). 


2.1 TAREA PLANA DE LA TEORÍA 
DE LA ELASTICIDAD 


Para la solución de una serie de tareas 


de la teoría de la elasticidad puede apli- 


carse el esquema de cálculo en forma de un 
estado tensional plano, o de un estado de- 
formacional plano, unificados bajo el nom- 
bre general de "tarea plana", la solución 
del cual nos lleva a la integración de la 
ecuación diferencial en derivadas parciales 


relativo a la función de tensiones: 


AA (cy) =0 (2.55) 


donde: 


A; Operador de Laplace en coordenadas 
rectangulares, 


Las tensiones normales y tangenciales 


7) oy Ty se expresan a través de la fun- 
ción p(x%,y) por las fórmulas: 


2 2 
0. op NE ; 
dy Y ay? 
T = Pe 
Ty — dxdy (2.56) 


si la ecuación (2.55) se resuelve para 
la región en la cual (- < y < %) se puede 
aplicar a ella una transformación de Fourier 


a través de la coordenada y representando: 


=P erre Ya 
Zn PO 


Fix, ¿) > 


encontramos en lugar de (2.55) la ecuación 


diferencial ordinaria: 
LE Ey rut) =0 
da* 


la solución de la cual: 


Pía) =U + Bayo +/% + (0 + Dx) d */* 
(2.57) 


Aquí A(E), B((E), C(E) y D (£) son fun- 
ciones arbitrarias de f determinadas a par- 


tir de las condiciones de borde de la tarea. 


Aplicando la transformación de Fourier a 


las igualdades (2.56) encontraremos las 


transformadas de las tensiones: 


0, (t) = 


(2.58) 


l 
were 
mw 
ty 
E 
de 


% (0) = 


PF) 
di? 


(2.59) 


> Se 3 o DE 
al A dy = 


(2.60) 


De esta forma, si es conocida la función 
de tensiones F(x,t) » las tensiones buscadas 
pueden ser determinadas a través de la in- 


versión de las expresiones (2.58) a (2.60) 


1 > a Ey 
Oy) ——f E *Fjtje Ya 
E e 


(2.61) 


09 2 A 
o (2 ,Y) = + === AIF (E)y Y y 


nan ». dy 
| (2.62) 
E (a dPE)  —iky 
7 y). a [7 LEE) ¿EY 
=y Var 0 
(2.63) 


En calidad de ejemplo concreto analice- 
mos el caso de aplicación de una carga ver- 
tical uniformemente distribuida q(y) a la 
superficie libre del semiplano. Orientan- 
do el eje (y) a lo largo del borde libre 
del semiplano y el eje x perpendicular a 
él en dirección hacia abajo, en la fronte- 
ra del cuerpo para x = 0 se deben satisfa- 


cer las siguientes condiciones de borde . 


ARS Ty =0 (2.64) 


Teniendo en cuenta que para la tarea 
analizada cuando x > las tensiones deben 
ser finitas, de (2.57) se obtiene que: 


Fly = Urbe (2.65) 


Transformando las condiciones de borde 


(2.64) tenemos: 


0 (0/5) = 0) 77, (0,8) =0 (2.66) 


donde a través de Q(f) se representa la 
transformada de Fourier de la función q (y) 
ar = —— jay ay (2.67) 


Sustituyendo en (2.66) las expresiones 
(2.58) y (2.60) obtenemos teniendo en cuen- 
ta (2.65) un sistema de dos ecuaciones al- 
gebraicas para la determinación de las 


funciones A(¿) y (E): 


- EA) = Q() 
El -/E/ A(E) + 5(£)] =0 


Determinando estos coeficien:es y susti- 
tuyendo los mismos en la ecuación (2.65) 
encontraremos la función P(x,£), y después, 
de acueráo con (2.61) a (2.63), las compo- 


nentes buscadas de las tensiones: 


== —— O bar 18/1002 EY 


x E ps 


(2.68) 
a 7 eo E =/E/x1E y 
o, TF La QU) 1 -/E/x]e dk 
(2.69) 
"A 00 7 Elxity 
ej = rd Q()e dé 
(2.70) 


Analicemos el caso particular de las 


Cargas en forma de fuerza concentrada apli- 


cada en el origen de coordenadas, esto es: 


qu) =P 5(y) 


Entonces, determinando de (2.67) Q(É), le 
encontraremos a partir de (2.68) y (2.70) 


los valores de las tensiones en el semies- 


pacio. 
3 
AS ES = 
(a + yy? 
2 
ls ZE 7 
y lo? +y) 


2.8 TRANSFORMACIONES INTEGRALES FINITAS 
DE FOURIER 


Como es conocido, bajo determinadas con- 
diciones la función arbitraria f(x) dada en 
el intervalo [0,a] se puede desarrollar en 
serie de Fourier por las funciones propias 
u (esk) de una cierta tarea de borde, o lo 
que es lo mismo, representar la función en 


forma de una serie: 


fu) = 3 FA, JubesA,) (2.71) 
donde : 

a 
FA) =f f (0) u(z,A Jdz (2.72) 


Se puede definir la fórmula (2.72) como 
la transfomación finita directa la cual 
transforma la función dada f(x) en una fun- 
ción FA), en ese caso, la fórmula (2.71) 


determina la transformación inversa. 


Si, por ejemplo, u (u,A) es una función 
de la ecuación elástica de la viga, esto es, 
una función dada por la ecuación diferen- 
cial: 


(2.73) 


y por los cuatro tipos de condiciones de 
borde: 


1. u(0) =u'"(0) =u(a) =u"(a) = 0 


2. u(o) =u'(0) = u(a) = u'(a) = 0 
3. u"(0) =u'"(0) = u"(a) = u'"(a) = 0 
4. u'(0) = u""(0) =u'(a) = u'"(a) = 0 
En el primer caso de condiciones de borde 


los valores propios serán: 


u (u/A,) = sen AZ ¿A == yama 


n 


2 


a 
F¿ A) da f(x) sen A Edo (2.74) 


fu) =,2, FSA) send,z 


y nosotros obtenemos un par de fórmulas, 
las cuales determinan la transformación-seno 


y finita de Fourier. 


El cuarto caso de condiciones de borde 


resuelve la siguiente ecuación: 


u (x,A,) = COS Ag (2.75) 
y entonces: 
Ñ 7 Aude 7 
E¿Aj) Ser! f(x) cos ye (2.76) 
a, $ 
fin) = a FA) cos A, 2 (2773 


Las fórmulas (2.76) y (2.77) determinan 
la transformación coseno y finita de Fou- 


rier. 


En el caso del segundo tipo de condicio- 
nes de borde tenemos que: 


TA) 
u(u,A ) = U(A,2) Sl TA vaz) 
y en el tercero: 

T(A_a) 


n 
u(x,A,) = S(,x) - TRA TA) 


Es evidente que la aplicación de las 
transformaciones integrales finitas de Fou- 
rier no da ninguna ventaja durante la solu- 
ción de la tarea en comparación con el mé- 
todo tradicional de desarrollo en series. 
Sin embargo, la misma permite realizar la 
discritización de la solución de las ecua- 


ciones diferenciales para regiones finitas, 


Le forma similar se introducen las 
transformaciones integrales dobles y fini- 
tas de Fourier, De esta forma, si la fun- 
ción f (x,y) está dada en la región que se 


cglcula[0<x <a, 0< y S bl], entonces: 


ab 
FAA? =/ Í f(x y) ue) u (y ¡A )dudy 


(2.78) 


Fic) = tb PALA) u(x,A,) U (y ¡A,,) 
(2.79) 


En el siguiente epígrafe será analizada 
la aplicación de estas transformaciones a 


la solución de tareas concretas. 


2.9 VIBRACIONES DE UNA PLACA RECTANGULAR 


Determinemos las formas de vibración de 
una placa rectangular, de rigidez cilin- 
drica D con dimensiones 4 y b, que se en- 
cuentra sometida a la acción de una carga 
dinámica arbitraria qlo yt). La ecuación 


de las vibraciones de la placa dada tiene 


la forma: 
2 
DAM fx y, t) om LEO y, t) 
dt* 
(2.80) 


En este caso Mes la masa de la unidad 
de área de la losa y A es el operador de 
Laplace en coordenadas rectangulares. 


En el borde de la placa pueden ser es- 
tablecidas diferentes condiciones de borde. 
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1. Apoyo articulado: 


a? 
w= 0 == -0, (r =0, x= =a) 
dx? 
0 
v=0 E O YESO by 
dy? 
2. Borde empotrado: 
du 
e Mr a (1 =0, x=) 
dw 
A DB) 
3. Borde libre: 
3% 0 “vn 
+ = ; E 
a? dy? da? 
y? 
+ (2 - y) — =0(0=0,%=a) 
dxdy? 
?w du 0 w 
Y =0; 
ay da? dy? 
3 
MA A 
dyda* 
4. Ligadura deslizable: 
3 
-=0 e (Tr =0, x= q) 
da? 
0w 0 
3 =0 =0 Y =D) Hi=%b) 
y dy? 
Para la solución de la ecuación (2.80) 


multipliquemos en ambos lados de la ecua- 
ción por el producto de las funciones 
u(x,An) UlyaA ds tenemos que u(<,A2) y 

u (y ,Am) son las funciones propias de la 
ecuación (2.73) bajo las correspondientes 
condiciones de borde en las ligaduras de la 


placa. Representemos: 


b 
f w(x yt) Uy (A) 7 (Y 1) dx dy 


o= pa 


Me NA 


(2.81) 


ab 
Or (5) =$ Gtrry rt) 4, (0/A,) Uy (Y A DÍA 
(2.82) 
Si integramos por parte obtenemos: 


3% (a y yt) a 
A Uy (5 A) Uy Y 1 9 da dy 


a b 
dE 
0.0 da 


¿e b 
= 10 S S w UY, drdy 


u! - 


A Po UM we 
o n xn xn 


= (ye 


LIL (2) 


a = xt 
do Ya Y A dy = A rm 


(2.83) 


Aquí se tuvo en cuenta que en correspon- 
dencia con las condiciones de borde la ex- 
presión dentro del corchete se anula para 


r=0,1=Q. De forma semejante: 


a b 4 
0 wir y yt) e ves, 
S/ 7 U,, (U eN) Y, Y ¿A )dxdy 


= AM y 


Att) 


(2.84) 


Integrando por parte, encontraremos: 


Diga 
S Ala y rt) u (e A) u(y,, )dudy = 


a 
f 
0.0 da dy? 


ab b u 
E n 9 Prat " 
f fu A dxdy + S [uo 4,4 Ue dy + 


a Nor Dn a 
+/ La une Y da + 


+ da, Enga Y? — a Lg 20: 
(2.85) 
La expresión obtenida muestra que en 
los casos en que el apoyo de la losa se 
corresponde con la primera y la cuarta 
forma de las condiciones de borde, en 
ellas se queda solamente el primer suman- 
do, y los restantes en estos casos se 
anula. Esto está vinculado al hecho de 


que en estas condiciones las funciones 
propias representan funciones senos y co- 
senos, por lo que obtenemos en lugar de 
(2.83) que: 


UD "ze 
7) ' 
LS Py rt) ulx,k ju(y,A )dedy = 
0.0 datay? n m 
2 2 
A q NE) (2.80) 


Teniendo en Cuenta las expresiones 
(2.31) hasta (2.86) la ecuación (2.380) se 
puede transformar en la ecuación diferen- 


cial ordinaria: 


vw. (t) 
Via ts y 


dt* 
(2.87) 


En este caso se representa: 


e 0 2 2 
a Al p As 3 a 


Por solución de la ecuación (2.37) te- 
nemos la función: 


Wan (t) = An cos Sm + Bi Sen Sp ca 


¡ya 
+ Wan E) 


La solución particular Wan CE) puede ser 
obtenida por aplicación directa de la 
transformada de Laplace a la ecuación 


(2.87). Como resultado obtenemos: 


ri q (pjeb* 


1 
Lic (t) = Si (2.88) 


2h 

¿oo 2 2 
LS po +S 
donde: 


Can (P) Transformada de Laplace de la 


función — Lam (5). 


Sin embargo, es posible no calcular 
Bram P? y en su lugar utilizando el teorema 


sobre la convolución, obtener directamente 


a partir de (2.88) que: 


1 t ; 4 
* = — K(t = rjiar 
We 50 70, (IKE 1) 
donde en su orden, si se utiliza la inte- 


gración a través del contorno, se obtie- 


ne que: ? 
1 Y +20 px 

k(0) =—_— S dp = 
21 . 2 2 

A EP Sam 
= sen £ zx 
mi 
Consecuentemente: 


É£ 4 

T (t) = Am cos Ba + Bón sen e 

1 t 
AN AS rjar 
rn nm 
US o 

Los coeficientes A y B se determinan de 
las condiciones iniciales de vibración de 
la placa. Si estas condiciones tienen la 


forma: 


9 Ub 
w(c,y,0) = fíie.y) y o 2 l ezo = g (1,4) 


entonces, transformando estos, obtenemos 
j isf dl 
que la función Mam E) debe satisfacer las 


condiciones iniciales. 
a 
Hg LO) o "4 Lezo P: en 


donde: 


ab 
F el S f(c,y) u(e,A,) u (y ¿A,,) dxdy 


ab 
=[[ 959) ule,,) uly,A,) deay 


Suponiendo ahora que las condiciones 
iniciales de vibración de la placa son 
nulas, entonces 2 = E = 0. Determine- 
mos ahora las vibraciones forzadas de la 


Placa en el caso de tener una fuerza con- 
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centrada p, aplicada en el punto EY Y 
que varía de acuerdo con la ley, 
Q(L,Yat) =P 8 (3 - y) Y y) 0.0%) 
(2.89) 
Entonces, sustituyendo (2.89) en (2.82) 
encontraremos: 
Can EY PAE UN y 1 0 1) 


y 


1 t 
as] am 7) sen an E = 1)de «Pa 


un 2 
pu A 


(1 = cos E) UA, ne, ) pa A) 


Consecuentemente a través de la fórmula 
12.79) obtendremos la expresión buscada de 
las formas de Vibraciones forzadas de una 
placa rectangular, sometida a la acción de 


una fuerza concentrada aplicada súbitamente, 


ESA (I=-cosS  t)u(A q.) 
ayer E A E a A) 
y abpD n=1 m=1 A 2 12 
7 m 


Esta solución corresponde al primero o 
cuarto caso de condiciones de borde. En el 
Primer caso, las funciones propias se deter- 
minan de acuerdo con (2.74) y en el segundo 


a (2.75). 


Si en el borde de la placa están dadas 
las condiciones de borde del segundo o del 
tercer tipo, entonces la solución de la 


tarea se complica considerablemente, 


Analicemos, por ejemplo, las vibraciones 
de una placa empotrada en el borde. En 
este caso, de la expresión (2.85) no se 


Puede obtener la expresión (2.86) y tenien- 


do en cuenta las condiciones de borde: 


eu 
da? dy? 


ab 
S f u,sA,) uy A, )dedy = 


ab ER 
.l S 5 Y dxdy = Pm 2) 
y la ecuación (2.80) se obtiene que: 


DOS, + A y (E) + 2 F,,, (8)] + 


dé A) 


+ pu =p 


dt? 


Para la solución de esta ecuación susti- 
tuyamos en la función Pi (6) la expresión w 
en forma de la serie. 


ola yrt) = EEN, (e A July A) (2. 90) 
nm 


Los coeficientes el se determinan por 
la fórmula (2.81), entonces: 


Pm E) = 22 A a) Arm (2.91) 


nm 


ab 
A el Í UA p5)u" (A, y Yu A, cu ¿Y )dxdy 


En resultado de la sustitución de la ex- 
presión (2.91) en (2.89) obtendremos un 
sistema infinito de ecuaciones diferencia- 
les ordinarias, resolviendo el cual deter- 


minamos los coeficientes Ma" 


4 a , 
na - O +22 rra + 


n m 
A E 
+ ——— =Q_ it) 
dt* sa 


Después de la determinación de los coe- 
ficientes JT la solución definitiva de la 
tarea la obtenemos a través de la fórmu- 


la (2.90). 


Capítulo 3 


TRANSFORMACIÓN DE HANKELL 


3.1 DEFINICIÓN DE LA TRANSFORMACIÓN 
DE HANKELL 


Durante la solución de ecuaciones dife- 
renciales en coordenadas polares en una 
región circular con cambio de la variable 
independiente r desde o hasta aplicamos 
la transformación de Hankell, la cual 
transforma la función dada f(r) en una 
función a, (A), de acuerdo con la fórmula: 


E, QM) El fir) rd, Arjdr  (0<A<oo) 


(3.1) 


para lo cual se supone que la función f (r) 


es tal, que se cumple que: 


57 Fo) VY dr <a 


En lugar de la igualdad (3.1) se utiliza 
también la nomenclatura 4, l f(r)l. La trans- 
formación inversa permite encontrar la fun- 
ción f(r) en el caso de ser conocida la 
función E, a). 

cn 
fr) = S E a) d,, (Ar)jadk (3.2) 

En las fórmulas presentadas dy (Ar) son 
las funciones de Bessell de orden 7 (y>-3). 
Para los casos axialsimétricos el indica- s 


dor Y = 0. 
3.2 PROPIEDADES DE IA TRANSFORMACIÓN 
DE HANKELL, 
3.2.1 Propiedad de linealidad a partir de 


A partir de (3.1) es evidente que: 
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£,Laf, (1) +Bf, (0)] =0x,[f (2)] + 


1 
£ iy 
+ Bx_1f, (01 
En este caso 4 y f son constantes y 


eE (r)] son las transformaciones de Han- 


kell de las funciones Pe (EJ: 142)5 


3.2.2 Propiedad de semejanza 


Determinemos la transformada de la fun- 


ción f(ar) donde 4 es un valor constante. 


A partir de (3.1) se ve que: 


fs Sis 
£,[F (07)] =/ flar)e d, A2rdr ds drag) 
(3.3) 


3.2.3 Derivación de la función original 


Determinemos la transformación de Han- 
kell de la función f*' (r). 


Durante la integración por partes se su- 
pone que el término fuera de la integral se 


anula, y se obtiene en este caso que: 
1%" (2)rd, Arjdr = Ls (Y - 
0 E 2Y 0 
a 109,_/ Ar) dr 


de esta forma: 


_—_— 


E 17 (r)] = E [(r - DE, 41 A) - 


MODE 02 (3.4) 


Determinemos la transformación de Han- 
kell de la función Af(r) donde A es el ope- 
rador obtenido del operador de Laplace en 
coordenadas polares, después del desarrollo 


de este a través de la coordenada de giro: 


a 


Integrando por partes, y suponiendo que 
se anula el término que está fuera de la 
integral, tenemos que: 


J Af(r)r d,, Arjdr .s fir, Ar]” - 


e , Vil O Y 
dy Ar) - Yo d, (Ar) hdr = “A E Ar) 
esto significa que: 
«LAf(0)] =- ”E, (Ar) (3.5) 
De esta fórmula resulta que: 


=. 24 
«1 ASP (1) A E, (Ar) (3.6) 


3.2.4 Transformacion de Hankell 
de la función r”' f (7) 


De la relación siguiente: 


dy 02) =Ard, (7) +, 001 (an 


y+1 


resulta que: 
-] á -1 
e Lr2fm) = ME, _,0) + E,,,0)] 21] 


3.3 DISTRIBUCIÓN DE LAS TENSIONES EN CUERPOS 
DE REVOLUCIÓN CON SIMETRÍA AXIAL 
La determinación de las tensiones en 
cuerpos de revolución que tienen un eje 
de simetría 2, puede ser transformada en 
la determinación de las llamadas funcio- 
nes de tensión y (r,2) de la ecuación biar- 
mónica: 
AA P(r,2) = 0 (3.7) 


Aquí Á es el operador de Laplace en 
coordenadas polares, el cual en consecuen- 
cia de la simetría y = 0. 


Si las funciones y (1,2) son conocidas, 
entonces las tensiones buscadas se hallan 
por las fórmulas: 


o. = FT IM 24 +2) 21, 


2 


0 
Ps A ] 
0, = . [AA pes + a (3.8) 
0 92 bs r dr ?*' A 
mb A a 2 y Ly] 
rz 7 A + JA p p 28 


donde : 


0 EG Tensiones en las direcciones ra- 
diales y circunferenciales. 


ez : Tensiones tangenciales. 


0. Tensiones en el plano perpendicu- 
lar al eje 2. 


A través de A y l se representan los 
coeficientes de Lamé, expresados a través 
del módulo de elasticidad E, y del coefi- 
ciente de Poisson Y por las fórmulas: 

E 


ii Ear 
1-21) 2(1 +7) 


Para la solución de la ecuación (3.7) 
apliquemos a ella la transformación de 


Hankell, representemos: 


96,2) =[7 vtr,2)10, Erdár 


Entónces en lugar de la ecuación (3.7) 
obtenemos la ecuación siguiente: 


d* 2,2 
— -¿*) pE,Z) =0 (3.9) 
dz 


la solución de la cual tiene la forma: 


+/ 2l El Z 


6 ,3) = (A + B2Z)e + (C + DZ)e 
(3.10) 


La función buscada de las tensiones 
p(r,z) se determina por la fórmula de in- 
versión: 

92) =/" 06,200, Ema 


(3.11) 


Sustituyendo (3.11) en la fórmula para 
las tensiones (3.8) representemos ellas a 


través de la función P(A,z2) por ejemplo: 


ES 3 
o,=fHara EL - 
0 da? 
E dl 

(A + 4) E” 10, Erdag (3.12) 
roer patas ames Ena 
PE AO dz? 1 

(3.13) 


Presenta también interés la determina- 
ción de los desplazamientos del cuerpo. 
Así los desplazamientos verticales se deter- 


minan por la fórmula: 


E AZ E 
0 dz? 4 


(3. 14) 


Las funciones arbitrarias que aparecen 
en las fórmulas (3.10)-(3.14) pueden ser 
encontradas de la satisfacción de las con- 
diciones de borde dadas en cada caso con- 


creto. 


Analicemos, por ejemplo, un cuerpo axial 
simétrico, limitado por el plano Z = 0, 


esto es, por un semiespacio. 


Representemos la intensidad de la carga, 
aplicada perpendicularmente a la superficie 
libre del semiespacio a través de q (r). En 
consecuencia de la limitación de la función 
de tensiones en el infinito, será necesario 
colocar en la expresión (3.10) los valores 
C=D=0. 
ra 2 = 0 se deben satisfacer las condi- 


ciones: 


Sobre la superficie de fronte- 


A A 


Su“tituyendo aquí (3.12) y (3..13) encon- 


Eramos: 


ABELLA + JE]? (E) = 


A (E) 


BE) = 2EIN2A + 0%], 


donde a través de Q(É) se representa la 
transformada de Hankell de la función q (r), 


esto es: 


q) = ¿que d, Ende 


y consecuentemente: 


9(,2) > 


(3.16) 


Analicemos una serie de casos particula- 
res de carga. Sea una fuerza unitaria con- 


Centrada actuando sobre el semiespacio, en- 


tonces: 
q (2) = (2n1)7* 5 (r) 
A partir de (3.15) se ve que: 


Q() = 1 


- 


y de (3.14), teniendo en cuenta (3.16) en- 
contramos, por ejemplo, el desplazamiento 
vertical de la superficie frontera del se- 


miespacio: 


y A + 2p Es 
u | ERA 1 d, (Er) de 


2=0 


a A + 24) 
A + pr 


De esta forma se Obtiene la conocida 
fórmula para la determinación de los des- 
plazamientos del semiespacio bajo la ac- 


ción de una fuerza concentrada. 


Si la carga está uniformemente distri- 
buida a través del área del círculo de 
radio r = a 
q(r) = q(r <a) 


Q) = qa Y, (ta) 


y de (3.14) encontraremos los desplazamien- 
tos de la superficie del semiespacio: 


rc id 


3=0 


3.4 VIBRACIONES LIBRES DE UNA PLACA 
INFINITA 


Analicemos las vibraciones libres d 
una placa lo suficientemente elástica, para 
que durante la determinación de los modos 
de vibración la placa se pueda considerar 
infinitamente grande. La ecuación de las 
vibraciones libres axialmente de la placa 
tiene la forma: 


wy,t) 


=0 


a? 1 a 2 
a Pa? wr, ,t) + 


at? 
(3.17) 
C?Dm*,me=ph, D = ER [12(1 - 1*)]”* 


donde: 
h: Espesor de la placa. 


RE/Y: Peso específico, módulo de elas- 
ticidad, coeficiente de Poisson 
del material de la placa. 


w(r,t): Forma de las vibraciones pro- 
pias: 


Supongamos que la condición inicial del 
movimiento de la placa es tal, que: 


w(r,0) me, Said =0 


t=0 


Aplicando a la ecuación (3.17) y a las 
condiciones (3.18) la transformación de 
Hankell con un Índice y = 0, obtenemos la 
ecuación diferencial ordinaria: 


CEWIESt) + A =0 
dt 


(391, 


con las condiciones: 


W(E,0) =€ , W, E 10) =0 (3.20) 


donde: 


W(E ,t) Transformada de Hankell de la 
función w(r,t). 


La ecuación (3.19) la solucionaremos 
con ayuda de la transformada de Laplace. 
Representemos: 


PE) SW ed, 


entonces a partir de (3.19) teniendo en 
cuenta las fórmulas (1.9) y las condicio- 


nes (3.20) encontraremos: 


Wip). ——— (3.21) 


p $ Es 


por la fórmula de inversión: 


ario 
eos. Hs 


aio poz 


Esta integral se puede calcular fácil- 
mente a través de una integral de contorno. 
Sin embargo, si se utilizan los ejemplos 
presentados en la fórmula (1.7) se puede 
directamente, a partir de la expresión 
(3,21), escribir que: 


WE ,t) =€ cos (ck*t) (3.22) 


La ecuación buscada de la forma de vi- 
braciones propias de la placa la encontra- 
remos por la fórmula de inversión para la 
transformación de Hankell: 


e(r,t) = ef” E cos (ch? 11d, (Enidé = 


E ———_— Sn 


2 ct E (t>0)r>0) 


3.5 TRANSFORMACIÓN FINIJA Di HANKELL 


Las transformaciones integrales y fini- 
tas de Hankell se aplican durante la reso- 
lución de ecuaciones diferenciales en coor- 
denadas polares para regiones acotadas o 
limitadas. Sea la función f(r) que existe 
en el intervalo finito (0 <r <a), la 
transformación finita de Hankell será la 


función: 
a 
PEL) =S fíejr pin E Pidr (3.23) 
donde: 
Es Raíz de la ecuación e, (2) =0 


La transformación inversa, en este caso, 
se escribe en forma de serie y determina la 


función f (7). 


J (E.r) 
pm = re) — (3.24) 
a E E 
y+1 91 


Si la ecuación se resuelve para una re- 
gión anular (b<r< a), la transformación 
directa y finita de Hankell toma la forma: 


a 
EE E Fra, (E¿mY, (ak) - 


- l, (E¿PI4a (az ¿Jar (3.25) 
donde: 
Y, (u): Funciones de Nieman de Índice y 
mientras la inversa: 
2:92 
Edy 60) a) 


Pr? > 2 2 
í e (E) > (E ¿b) 


(9, EL0%, 60) =D, (2%, (E 2)] 
En tareas prácticas, frecuentemente se 


exige determinar la transformación finita 
de Hankell de la función Af (r) 


donde: 


A: Operador a Laplace en coordenadas po- 


lares: 


En el caso de una región que sea circu- 


lar (0 <Y <a), integrando por partes, ob- 


tenemos: 

a 

S Afimr Y, (E ride = Pla) dla) — EFE) 
(3.26) 

donde: 


FEs): Se determina por la fóxrmu- 


la (3.23). 


Para una región anular (hb <r< a): 


a 
[ mad, E. , Ea) - 
d, (Ea) 
MODA de 
E d+ E¿DY, (E r)Jar =-¿ 7, 72) F(b) + 
2 2 E 
FO) - EFE) (3.27) 
donde: 


F,s Se determina por la fórmula (3.25). 


3.6 DISIKIBUCIÓN LEL CAIOR EN UN CILINDRO 
DE LONSITUD INFINTIA 
Determinemos la distribución de la tem- 
peratura en un cilindro circular Maazo 
(0 < r <€ 0) de gran longitud, si su super- 
ficie se somete a una función de tempera- 
tura, variable en el tiempo f(t). Orien- 
temos el eje z a lo largo del eje del ci- 
lindro y supongamos que la temperatura 
no varía a lo largo del mismo y que se en- 
cuentra simétricamente distribuido con re- 
lación a él. Supongamos, además, que la 
distribución inicial de la temperatura se 


determina por la función f, (r). 


La terperatura T(Y,t) en la sección del 
cilindro a una distancia r del eje y en un 
momento de tiempo + se determina a partir 


de la solución de la ecuación en derivadas 


parciales: 
1. APM) ¡ET 1 9T(P,t) 
K t SÍ dp? $ r dr ) 


(3.28) 


donde : 


k: Coeficiente de transmisión del calor 
del material del cilindro. 


Las condiciones iniciales y de borde de 


la tarea tienen la siguiente foma: 
T(2,t) = F(t) (r =a) t>0) y T(r,0) = 
s. La (2) 


Teniendo en cuenta que la distribución 
de la temperatura es axial simétrica, en- 
tonces, aplicando a la ecuación (3.28) la 
transformación de Hankell en el caso en 


que Y = 0 y considerando (3.26) obtenemos: 


drE,t) : 
ió KEz FE jet) ka EY, E¿2)f(t) 


(3.29) 
La condición inicial transformada: 
a 
Ei Ez) +] rf, (md, (E,P)dr (3.30) 


La solución de la ecuación (3.29) bajo 
la condición inicial: 


FEL¿0) = E, E) 

se escribe en la forma: 
=kg 1? 

FE,¿% =F le + 


t KE, (t) 
+ ka Ed, EDS fítr)e dí 


Por la fórmula de inversión (3.24) ob- 
tenemos la función buscada de distribución 
de las temperaturas siguientes: 


da tr, 2.) 


Pr t) 2 4... 8) E 
a PAN: 
E az 


Si, por ejemplo en el momento inicial de 
tiempo la temperatura se iguala a Cero, 
esto es f, (P) = 0 y F, (E;) = 0, entonces: 


. EME E =kE? (t) 

2: sy EU 0 27 > 1 P 

10 Er lim de 
13.31) 


Sustituyendo en (3.31) la distribución 
dada de la temperatura sobre la superficie 
del cilindro encontraremos su valor en el 
interior del cilindro en un cualquier ins- 


tante de tiempo lt. 
Por ejemplo, si sobre la superficie del 


cilindro la temperatura es constante: 


fit) =f 


entonces: 


2 
E (E,1) kE, t 
¿, (Ea) 
(3.32) 


La fórmula 3.32 da la distribución de 
temperatura en el interior del cilindro, 


si sobre su superficie ella tiene un valor 


Capítulo 4 


TRANSFORMACIÓN DE MELLIN 


4.1 DEFINICIÓN Dk LA TRANSFORMACIÓN 

DE MELLIN 

Para la solución de ecuaciones diferen- 
ciales en coordenadas polares, dadas en 
una región en forma de cuíia (0 < yr < 00), 
(a <0 <a) aplicamos la transformación 
integral de Mellin la cual se determina 
por el siguiente par de fórmulas. Se llama 
transformada de Mellin de la función fir) 
a la función: 


F(s) = 3 E 


La transformación inversa se determina 
por la fórmula: 

PE 

4 + 


—_— f 


zz 


s+1 


Físjr "de 


0-00 


fí(r),= (4.2) 


Aquí tal y como en (1.4), la integral 
se toma a lo largo de la recta Res = 0 > Sor 
que garantiza su convergencia y representa 
el límite de la integral a lo largo del 


segmento 0 - 1b, 0 +1b, parab>o, 
4.2 PROPIEDADES DE LA TRANSFORMACIÓN 

DE MELLIN 
4.2.1 Propiedad de linealidad 

De la fórmula (4.1) se tiene que si 
E (r) y E (r) son dos funciones que poseen 
las transformadas EME Y F, (8), entonces: 
Flaf, (1) + Bf, (2)] = ar, (s) + BF, (s) 


(4. 3) 
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4.2.2 Propiedades ae semejanza 


Para la función f (ar) tenemos que: 


co - 1 des s-1 £(s) 

SF are? da = TE Fu du 
[04 Q 

(4.4) 


donde: 


F(s): Se determina por la fórmula (4.1). 


4.2.3 Transfomación de la función rf (rm) 


La transformación de Mellin de la fun- 
ción rf (2) será la función F(s + 4). 
200 (3) 
Evidentemente / r fire “dr=PF(s +0) 


(4.5) 


4.2.4 Derivación de la función original 


Determinemos la transformación de Mellin 
de la función pa (r). Paran= 1 te- 
nemos: 


Sao = fr)re72 JE 


0 
-(s - 2) Sr ar = -(8 - 2)F(8s - 1) 


Se supuso que la función f(r) es tal, 


que el miembro fuera de la integral dentro 


de los límites de integración señalados se 
anula. De forma semejante, integrando por 
partes se obtiene: 

IS 8-2 n ys 
S fm dr = (1) P(e-n) 1, (6-1=1) 


(4.6) 


Si la función dada f nm) (7), entonces: 


n 
(ur (s) ¿H, (r+s-¿=1) 


md 
(4.7) 


4.3 TAREA PLANA DE LA TEORÍA 
DE LA.ELASTICIDAD PARA UNA REGIÓN 
EN FORMA DE CUNA 


Utilicemos los resultados obtenidos an- 


teriormente durante la tarea de determina- 


ción de las tensiones en una cuña plana 
(OGÍT<, « <0 <a) sometida a una 
carga uniformemente distribuida (q) perpen- 
dicular a ambas fronteras de la cuña 

(fig. 4.1) 


Figura 4.1 


Las tensiones interiores que actúan en 
el interior de la cuña en las direcciones 
radial y anular O, Y % Se expresan a tra- 
vés de la función de tensiones Y por las 


fórmulas: 
2 
as = E, dy De e (4.8) 
5 rr dr 
2 
0) = 2 (4.9) 
ir 
y las tensiones tangenciales: 
e 0 1 dp 
e E E (5570) 


en su orden, las funciones de tensión pue- 
den ser determinadas de la solución de la 
ecuación: 


bMp=0 (4.11) 


donde: 


A: Operador de Laplace el cual en coor- 
denada polares tiene la forma: 


ee a? E e a? 
dy” E $ y 29? 


En las fronteras de la cuña cuando 
0 = ta deben satisfacerse las siguientes 
condiciones: 
q(0<«r <a) 
EG m0 y, (4.12) 
0(r > 0) 


Escribamos la ecuación (4.11) en forma 


desarrollada: 
4 3 2 
LD 
in dy? r”r ¿e 
1 a 4 a? 2 a? 
+ y — _—— - —— + 
,» 292 r”  dyo0? 
2 at a* 
+ + — lp (7,0) = 0 
y? ar*a0? y 20 


y para su solución aplicaremos la transfor- 


mación de Mellin. 


Con este objetivo multiplicando ambas 
s+2 A 
partes de la ecuación (4.13) por Y e in- 
tegramos a través de r desde cero hasta ”, 


Representando: 


Ps 0) le pr 031 ?dy (4.14) 


Sobre la base de la fórmula 4.7 obte- 


nemos: 


e 3%ty pa+2 


S ; dr = (s-1)s (s8+1) (s+2) P(s 4) 
dr 

oo y s+2 

S o el dr = -(s-1)s(8+2) P(s 0) 

or 

oo 2 

A IC e (6d A 

a? 

¡— Z y "2 = -(8-1) P(s/0) 


Utilizando las expresiones obtenidas, 


tendremos en lugar de la ecuación (4.13) la 


ión dif > ' Pe ” 
ecuación diferencial ordinaria a Al A $ (8-10 
4) = LL (S+1) sen (s+1)a cos (s-1)0 = 


2 5 
dls 0) A E Pots 0) 4 sis” - 1)(s sen 2 + sen 254) 
de? dg? : 
(8-1) sen (s-1)a cos(s+1)9] 
TAS TE A A Cd 10 
AC 1) (s 2% + sen 250) 
+ (8 - 1 b(8 0) =0 EY > (4.16) 


La solución de la anterior es la función: Ahora a través de la fórmula (4.2) de- 


terminamos la función de tensiones: 


P(s,0) = A senís - 1)0 + 8 costs — 1) 0 + 
4 O+To0 
p(r,0) ES 
oo 


, =5+1 
+ C senís + 1) 0 +D cosís + 1)0 (4,15) ds (4.17) 


Las constantes arbitrarias se determinan E á 
Calculemos las tensiones tangencia- 


a partir de las condiciones de borde en las as ; 
E les Trg" “on este objetivo, sustituyendo 


í de 1 na, 1 les t ié 
ronteras de la cuña, las cuales también es UT] en (4.10) obtenems: 


necesario de acuerdo con la transformación 


de Mellin resar a través de la fun- ario 241 
6 E, n= Ip  _£ sen(er)a sen(s-1)0 
ción P(s,0). Prestemos atención al hecho AA 
q de si ae O ¿90 An S sen 2% + sen 220) 


que como consecuencia de la simetría de las 


cargas los coeficientes A =C = 0. Para - senís-1)a_sen(e!1)8] Mo 


la determinación de los restantes coeficien- pet (8 sen 2% + sen 280) 
tes, sustituyamos (4.9) y (4.10) en la ecua- 


ción (4.12) y obtenemos que parad =1tú: (4. 18) 


Durante el cálculo de la integral de 
q (0<rc< aq) > 
0% contorno en la ecuación (4.18) destaquenos 


li) 
+ LH Oy 


ap? el hecho de que la función bajo el signo 
E O (a <p) 
de integral tiene un polo en el punto 


s=0. 


Multiplicando ambas partes de las lggal= Consecuentemente escogemos el contorno 


E SS z 
dades obtenidas por 7”, integremos la misma de integración, compuesto de una recta que 


desde cero hasta infinito e nEBO, ER va a lo largo del eje ficticio con la ex- 


cuenta la ecuación (4.7) obtenemos las si- cepción del punto singular s = 0 y la cir- 
guientes ecuaciones para la determinación cunferencia de radio R es el semiplano de 

j j D. ; 
de las constantes arbitrarias By la izquierda (fig. 4.2). 


B cosís - 1) a + D cosís + 1) a = tu 


-1 
=- [s (s? -= 1)] ee Sr 


B(ls - l)senís — 1)a + y 
+ D(S + 1)sen(s + 1)Ja =0 


Determinando B y D y sustituyendo sus 
valores en la ecuación (4.15) obtenemos : Figura 4.2 
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Realizando los cálculos correspondien- donde: 
tes encontraremos que: 


Wís, 0) =[%0 (o, 0) 1? gn 


Ps En f”. sen (4-9 ) senh (440 )u- A 
r0 Tr Y u sen 2% + senh 2ua la solución de la ecuación (4.21) tiene la 
forma: 


. —sen (4+8) senh (4-9 )u cos (u 1n Edu 


u sen 2% + senh 244 W(s, 0)= A sen (s - 1)0 + B cosís - 1) 0 + 


* 
De forma semejante se determinan tam- + C sen (5 + 1)0 +Dcos(s + 1)0 +% (5,0) 


bién las otras tensiones. Para ello sus- ic 
tituimos (4.17) en (4.8) y (4.9) y reali- donde: 

zando luego la integración a través del W*(s, 0): Solución particular de la 
contorno al igual que en el caso analizac ecuación (4.21) y los restan- 
se obtienen las expresiones buscadas en tes sumandos representan la 
forma de integrales sobre el eje real. solución de la ecuación homo=- 


genea. 
4,4 FLEXIÓN DE LA LOSA EN FORMA DE CUNA 


Para la obtención de la solución parti- 
En calidad ae segundo ejemplo determine- 
cular de la ecuación (4.21) le aplicaremos 
mos las flechas de una placa en forma de 


a ella la transformación de Laplace. Si 
cuña w(r,0) («a <0 <a, O<r <*%) someti- 


, representamos: 
da a la acción de una carga arbitraria 


q(”, 0) perpendicular a su plano (fig. 4.3). 5 co 
; W (s ,p) ed Wi 0 go 
Q(s Pp) = ie Qí.0)e Pg 


obtenemos de (4.21) teniendo en cuenta la 
fórmula (1.9): 


A 
La ecuación de la flecha de la losa en A y 


coordenadas polares tiene la foma: ; = 
Las funciones W(sS,p) es el producto de 


DAA 6 (1,0) = g (1,0) (4.19) la transformada de Laplace de dos funciones, 
consecuentemente y de acuerdo con el teore- 


si representamos: 
ma de la convolución la solución particular 


l 8-2 buscada: 
Q(s, 0) a q (2, 0)” dr (4.20) 
* on 
entonces, tomando en cuenta los resultados W (80) dl Q(s,T) k(s,0 - r)dT (4.25) 
del ejemplo anterior obtenemos en lugar de 
donde: 

(4.19) la ecuación diferencial ordinaria: 
A Pp k(8,2%) . a Le dp 

Wís, 0 dW(s, 0 3%) 5 GT E 
E a y A E quis pp + pr? 

d0 do? 

+ (8% - 1)7W(8, 0) =Q (s, 0) (4.21) (4.281 
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Las raíces del denominador de la función 
bajo el signo de integral en la ecuación 
(4.24) son: 


Pp e DS = il (p + 1) 


1,2 3,4 


Por eso, realizando la integración a 


través del contorno, obtenemos: 


sen (3+1)x ] 


k(s,2) = E Inez z a 


ed: , 
(4.25) 

De esta forma, las fórmulas 4.20 y 4.25 
permiten determinar la solución particular 
(4.23). 
centrada p aplicada en el punto de la losa 


Sea, por ejemplo, la fuerza con- 
de coordenadas = 0, , en este caso: 


q(r/0) =prs (» - r)) 5(0 -6,) (4.26) 


donde : 
6 (2): Función delta. 


Sustituyendo (4.26) en (4.20) obtenemos: 


Q(8,0) = PD? y 


8+1 
o. 50 -8,) 


Consecuentemente, de acuerdo con la fór- 
mula (4.23) y para0 <0, 


E 0 
2 l sen s(06 -0,) s 
2Ds(8g - 1) 
Wis/0) = <» cos (0 - 0,) - se cos 80 -0,), 
+ sen(0 - 0,)] 
(4.27) 


Las funciones arbitrarias en (4.22) se 


— => 


. determinan a partir de las condiciones de 
borde. Supongamos que la losa se encuentra 
empotrada rígidamente a lo largo de los 


bordes 4 = 0,0 =0 se obtienen: 


m(r,0) = uy (7,0) = w(r,a) = “y (r,a) 
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Las condiciones de borde transformadas 


tienen la forma: 


W(s,0) = hy (Ss ,0) = W(s,0) = W¿(s,2) (4.28) 


Sustituyendo (4.22) y (4.27) en (4. 28) 
obtenemos un sistema de cuatro ecuaciones 
para la determinación de las funciones ar- 


bitrarias buscadas, la solución del cual da: 


(s+1)f, (a,s) f, (4,8) 


A(s) = “ms 20) - Soo 
donde : 
0O(8) = sendas - s* sendas; 

f, (a,s) = W (a) (s cos 54 sena + 

+ sen Q cos 04) - Wo (4) sen sa sena; 
f, (0,5) = W¿(a) (s cos sa sena - 

- sen 34 cos 4) - 
- (s? 1)” (a) sen sa sen Q; 

también: 


w(a) y dy (4) se determina por la fórmula 


(4.27) cuando 0 =a, 


Ahora, de acuerdo con (4.22) la función 
w(s,0) es conocida y para encontrar la fun- 
ción verdadera de la flecha es necesario 


invertir la expresión (4.22): 


o+to SS 
ws jr “ds (-o, <R¿0<0, ) 


Escogiendo por la vía de integración el 


eje imaginario 0 = 0, obtenemos: 


AN CTE DES 


zi 
> senta s-s* sen a 


sen 8 Ócos 0)f, (a,s) + 
ic 


senta s-s*sen a 


+ sen 80 senó f, (0,8), 


€Q-———_ o a. Y) 


senta s-8* sena 


“Eg 


(4.29) 


Conociendo la función de las flechas, 
es fácil encontrar las solicitaciones de 
cálculo en la losa. 


Así por ejemplo, determinemos el mo- 
mento flector en la entikación de la 
losa cuando 0 = 0: 


a partir de (4.29) obtenemos: 


1-8 
Mtro) = - Al a pi O 
n *- cie 0(a,5) 


y teniendo en cuenta (4.27): 


M(r,0) = 


Pr, +leo sen 0, sen 8% sen s(0-6,) - 


= - —— Ká___—_—_—_——_—_—_——.- 
2ntp”., 2 2 2 
—Loo sen 84 - 8 sen Q 
- 8 sen senia4,) sen sí o 
. 2 (G> de 
sen 8 - 8 sen* a 


Pasando a las variables reales, obtene- 
mos la siguiente expresión para el momento 
flector en una losa empotrada en forma de 


cuña cargada con una fuerza concentrada: 


BL e 
M(r,0) | [ sen 0, senh ua sen u (ab, )- 


- u sen A sen(a-0,) senh u 0,]. 4 


r 
cos (u ln 


e 


— du 
sen? ya - 14 seña 
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